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ABSTRACT 


The  minimum  principle  for  velocity  fields  and  the 
maximum  principle  for  stress  fields  for  a  Bingham  solid  are 
extended  to  more  applicable  forms  to  obtain  upper  and  lower 
bounds  for  the  rate  of  work  done  by  the  unknown  surface  tract¬ 
ions.  These  extensions  are  applied  to  the  problem  of  pure 
torsion  of  a  square  bar  and  close  bounds  are  obtained  for  the 
torque  over  a  range  of  material  properties.  In  addition, 
other  approximate  solutions  are  presented.  The  torques  ob¬ 
tained  from  these  solutions  are  compared  with  the  upper  and 
lower  bound  torques  and  are  found  to  be  in  satisfactory  agree¬ 
ment  with  the  bounds . 

An  experimental  investigation  of  the  pure  torsion  of 
circular  cylinders  of  commercially  pure  lead  is  presented. 

The  results  obtained  show  that  commercially  pure  lead  behaves 
according  to  the  Bingham  visco-plastic  model  for  a  certain 
range  of  rate  of  twist  per  unit  length,  and  beyond  a  critical 
region,  the  torque  is  independent  of  rate  of  twist  per  unit 
length  up  to  the  maximum  rate  used  in  the  experiments.  The 
maximum  strain  rate  produced  in  the  experiments  was  0.36  per 
sec.  A  brief  investigation  of  the  shortening  of  pure  lead 
circular  cylinders  during  torsion  is  also  presented  and  the 
results  are  discussed. 
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NOMENCLATURE 

a 

=  semi-width  of  square  cross-section  in  the  torsion 
problem 

ei  j 

=  strain  rate  tensor 

0  0 
xz  ’  yz 

=  shear  strain  rate  in  the  torsion  problem  (*)+ 
e  e 

xz  yz  j .  .  ,  .  ,  , 

e  ,  e 

xz 9  yz 

.  ,  -r1 —  dimensionless  shear  strain  rate  m  the 

0a  9a 

torsion  problem 

F  . 

l 

=  body  forces  per  unit  volume 

I 

=  J  2e^1e^1  a  strain  rate  invariant  (*) 

I 

=  .  dimensionless  strain  rate  invariant  in  the 

0a 

torsion  problem 

J 

=  J  a  stress  invariant 

k 

=  yield  stress  in  shear 

l 

=  length  of  the  bar  under  torsion  (*) 

l 

=  —  dimensionless  length  of  the  bar  under  torsion 

nx>  ny 

=  x  and  y  components  of  the  normal  vector 

Si  j 

=  deviatoric  stress  tensor 

s  ,  s 
xz  *  yz 

=  shear  stresses  in  the  torsion  problem  (*) 
s  s 

xz  yz  .  ,  ,  .  . 

xz *  yz 

,  ,  —f —  dimensionless  shear  stresses  in  the 

k  *  k 

torsion  problem 

S 

=  surface  of  volume  V  (*) 

S 

g 

=  ~2  dimensionless  surface  in  the  torsion  problem 
a 

ST 

=  the  portion  of  S  on  which  surface  tractions  are 
prescribed 

+ 


An  asterisk  at  the  end  of  a  description  of  symbol  means  that, 
in  Sections  4.4  to  4.7,  the  symbol  is  used  in  place  of  its 
barred  symbol  unless  specially  noted. 
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the  portion  of  S  on  which  velocity  is  prescribed 

torque  in  the  torsion  problem  (*) 

T 

-tt—  dimensionless  torque  in  the  torsion  problem 
a  k 

components  of  a  stress  vector  applied  on  a  surface 
velocity  components  in  the  x  and  y  directions 
respectively 

stress  function  in  the  torsion  problem  (*) 

U 

—  dimensionless  stress  function  in  the  torsion 
problem 

components  of  a  velocity  vector 

volume 

V 

—  dimensionless  volume  in  the  torsion  problem 
a 

warping  velocity  in  the  torsion  problem  or  velocity 
component  in  the  z  direction  (*) 
w 

- — ~  dimensionless  warping  velocity  in  the  torsion 

0a 

problem 

coordinates  in  a  rectangular  coordinate  system  (*) 

X  v  z 

— ,  —  dimensionless  coordinates  in  the  torsion 

a’  a*  a 

problem 

yield  stress 

» 

a  dimensionless  parameter  in  the  torsion  problem 
parameter  used  in  the  Ritz’s  method 
coefficient  of  viscosity 

coefficient  of  viscosity  in  the  linearized  constitutive 
relationship 
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=  angle  of  twist  in  the  torsion  problem 
=  rate  of  twist  in  the  torsion  problem 
=  modulus  of  rigidity 
=  stress  tensor 
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s  +  s  magnitude  of  the  shear  stress  vector 
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in  the  torsion  problem 
=  speed  of  twist  in  rpm 


NOTE:  For  the  equations  which  are  in  tensor  notations,  the 

summation  convention  is  implied. 
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CHAPTER  I 


INTRODUCTION 


Many  solids  do  not  behave  elastically  under  load. 

Even  if  a  solid  is  elastic,  it  is  so  only  under  a  certain  range 

of  load.  The  deviation  from  elastic  behavior  in  a  solid  is 

generally  termed  "inelasticity" .  Inelasticity  in  a  solid  is 

caused  by  dissipation  of  mechanical  energy  by  the  solid  in 

three  ways.  These  three  kinds  of  inelasticity  are:  (a)  an- 

elastic  deformation  -  caused  by  the  dissipation  of  mechanical 

energy  through  the  interaction  with  flux  of  non-mechanical 

energy  such  as  thermal,  electrical,  and  magnetic  energy; 

(b)  viscous  response  -  caused  by  the  dissipation  of  mechanical 

energy  through  the  viscous  interaction  between  elements  of  the 

solid;  (c)  general  plastic  response  -  caused  by  the  dissipation 

* 

of  mechanical  energy  through  solid  friction  due  to  slip.  (1) 

In  a  "real"  solid,  all  three  forms  of  inelasticity 
are  present  to  some  extent  under  various  conditions.  This 
thesis  is  concerned  with  a  solid  for  which  the  second  and  third 
kinds  of  inelasticity  mentioned  above  dominate.  This  solid 
is  known  as  a  "visco-plastic"  solid  or  a  Bingham  solid. 

In  processes  such  as  hot  extrusion,  hot  forging,  and 

* 

Numbers  without  decimal  enclosed  by  brackets  designate  the 
references  listed  in  the  Bibliography. 
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hot  rolling,  the  strain  rate  dependent  properties  of  the  metal 
usually  influence  the  load  required  to  produce  deformation. 

The  extrusion  of  pure  lead  at  room  temperature  is  an  example 
of  such  a  process.  One  method  of  considering  the  rate  effects 
is  to  use  the  Bingham  solid  model  to  represent  the  metal  and 
the  purpose  of  the  experimental  part  of  this  thesis  is  to 
investigate  how  closely  the  behavior  of  commercially  pure  lead 
under  torsion  at  room  temperature  approximates  to  that  of  the 
ideal  Bingham  solid.  Preliminary  research  studies  of  the  rate 
effects  are  usually  done  with  pure  lead,  because  the  apparatus 
for  performing  experiments  with  hot  metal  is  quite  elaborate. 

In  the  theoretical  part  of  this  thesis,  extensions 
of  the  extremum  principles  for  a  Bingham  solid,  some  approxi¬ 
mate  methods  of  solution  of  the  problem  of  pure  torsion  of  a 
visco-plastic  prismatic  bar  with  a  simply-connected  cross 
section,  and  some  general  theories  of  plasticity  derived  for 
the  Bingham  solid  are  presented. 

It  is  hoped  that  the  investigation  of  some  methods 
of  analysis  for  the  problem  of  pure  torsion  will  give  some  in¬ 
dications  of  possible  methods  of  analysis  for  other  visco¬ 
plastic  boundary  value  problems,  and  the  experimental  investi¬ 
gation  of  pure  torsion  will  lead  to  a  possible  method  of 
determining  viscosity  coefficients  of  visco-plastic  materials. 
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CHAPTER  II 


GENERAL  THEORY  OF  VISCO-PLASTICITY 


2.1  Brief  Historical  Background 

While  the  development  of  the  theory  of  plasticity 
has  progressed  rapidly  during  the  last  twenty  years,  the 
interest  in  visco-plasticity  has  grown  only  in  the  past  few 
years.  In  earlier  years,  visco-plasticity  was  of  interest 
mainly  to  rheologists,  but  in  recent  years,  visco-plasticity 
has  also  been  developed  as  a  branch  of  applied  mechanics  that 
is  concerned  with  a  solid  that  possesses  viscous  properties 
when  in  the  plastic  state. 

As  early  as  1922,  a  model  for  a  visco-plastic  solid 
subjected  to  simple  shear  was  proposed  by  Bingham  (2) .  For 
this  model  the  shear  stress  due  to  a  certain  shear  strain  rate 
is  the  sum  of  the  shear  stress  for  a  Newtonian  fluid  and  the 
shear  stress  for  a  plastic  solid  due  to  the  same  shear  strain 
rate.  A  solid  that  behaves,  in  simple  shear,  like  Bingham 's 
model  is  now  known  as  a  Bingham  solid.  In  1934,  Grimshaw  (3) 
attempted  to  consider  the  rate  effect  on  the  plastic  deforma¬ 
tion  of  a  right  circular  cylinder  of  soft  metal  under  axial 
compression.  Grimshaw  assumed  the  stress  to  be  the  sum  of 
the  stress  in  a  Newtonian  fluid,  proportional  to  strain  rate, 
and  the  stress  from  an  elastic-type  relation,  proportional  to 
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strain.  The  modern  theory  of  plasticity  shows  that  the  use 
of  an  elastic-type  or  deformation  type  of  relation  for  plastic 
deformation  is  not,  in  general,  correct. 

Apart  from  Grimshaw's  work,  practically  all  other 
papers  on  visco-plasticity  have  been  based  on  the  generaliza¬ 
tion  of  the  model  proposed  by  Bingham.  A  constitutive  relation 
for  the  Bingham  solid  was  formulated  by  Hohenemser  and  Prager 
in  1932  (4) .  The  Mises  yield  criterion  is  used  in  this  model; 
consequently,  the  constitutive  relation  is  non-linear.  Prager 
in  1961  (5)  presented  a  theory  of  visco-plasticity  based  on  a 
piecewise  linear  yield  criterion. 

In  1946  and  1947,  Oldroyd  (6)  formulated  equations 
for  the  plastic  flow  of  a  Bingham  solid  and  solved  some  two- 
dimensional  boundary  layer  and  rectilinear  flow  problems, 
though  mainly  from  the  rheological  point  of  view. 

The  first  general  results  for  a  visco-plastic  solid 
were  obtained  by  Prager  in  1954  (7) .  Prager  derived  extremum 
principles  analogous  to  the  principles  of  minimum  potential 
energy  and  minimum  complementary  energy  in  classical  elasticity 
and  plasticity.  Extensions  of  these  extremum  principles  to 
more  applicable  forms  are  presented  later  in  this  thesis. 

Very  few  problems  have  been  solved  for  the  visco¬ 
plastic  model.  Carlson  (8)  considered  the  axial  compression 
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of  a  disk  without  accounting  for  inertia  effect  and  using  a 
piece-wise  linear  yield  criterion.  Haddow  (9)  solved  the 
same  problem  with  inertia  effect  considered  and  using  the 
Mises  yield  criterion.  Appleby  and  Prager  (10)  applied  the 
linearized  constitutive  relation  to  solve  the  flow  problem  of 
a  thin  circular  plate  subjected  to  a  uniformly  distributed 
transverse  load.  Appleby  (11)  also  considered  the  visco¬ 
plastic  expansion  of  a  long  thick  walled  tube  subjected  to 
internal  pressure. 

Ziegler  (12)  proposed  four  models  for  dissipative 
media.  Two  of  the  models  may  be  considered  as  visco-plastic, 
one  a  visco-plastic  solid  and  the  other  a  visco-plastic  liquid 
according  to  the  classical  definitions  of  a  solid  and  a  liquid. 
Ziegler  solved  a  simple  flow  problem  for  the  dissipative  media 
and  obtained  two  limit  theorems. 

2,2  The  Bingham  Solid 

The  Bingham  solid  is  so  far  the  most  common  visco¬ 
plastic  material  that  researchers  have  investigated.  In  the 
modern  literature,  the  term  Bingham  solid  refers  to  a  solid 
that  behaves  as  follows:  When  the  stress  state  does  not  violate 
the  yield  criterion,  the  solid  is  rigid,  that  is,  the  strain 
and  strain  rate  remain  zero.  When  the  state 
of  stress  exceeds  the  yield  stress,  visco-plastic 
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flow  with  no  workhardening  occurs.  Then  the  stress  intensity 
is  the  sum  of  two  parts:  the  stress  obtained  from  the  same 
strain  rate  applied  to  a  rigid  perfectly  plastic  solid  and 
the  stress  obtained  from  the  same  strain  rate  applied  to  an 
incompressible  Newtonian  fluid.  A  mechanical  model  of  this 
solid  is  shown  in  Fig.  2.1. 


Viscous 

Damper 


Fig.  2.1  Mechanical  Model  for  Bingham  Solid 


Mathematically,  the  stress  components  are 
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where  f(s. .)  is  the  yield  function,  and  (s. .)  and  (s. .)  are 
1  1  p  D  ^ 

the  stress  deviators  for  a  perfectly  plastic  solid  and  a  Newtonian 
fluid,  respectively,  subjected  to  the  given  strain  rate  . 

Using  the  Levy-Mises  relation  for  rigid  perfectly  plastic  solids 
(13) ,  in  which  the  Mises  yield  criterion  is  used, 
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Therefore,  the  constitutive  relation  for  Bingham  solid  is 


s^j  Indeterminate  if  1  =  0  that  is  if  J ^  k  (2.5) 
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if  I  f  0  that  is  if  J  ^>k  (2.6) 


where  J  =  s .  .  s .  .  is  the  square-root  of  the  second  invar- 
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iant  of  the  deviatoric  stress  tensor  and  I  =  s/2ev ^^i  is 
the  square-root  of  four  times  the  second  invariant  of  the 
deviatoric  strain  rate  tensor.  To  invert  relation  (2.6), 
both  sides  are  squared.  Thus, 


J  = 


_j_ 
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(2.7) 


then, 


I  = 


J  -  k 


(2.8) 


After  substituting  (2.8)  in  (2.6)  and  rearranging,  the  inverted 
constitutive  relation  is 
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2.3  Bingham  Solid  and  Rigid  Perfectly  Plastic  Solid 

In  this  section,  the  major  difference  between  the 
Bingham  Solid  and  the  rigid  plastic  solid  is  discussed.  A 
rigid  plastic  solid  exhibits  no  rate  effects;  the  current  state 
of  stress  defines  the  ratios  of  the  components  of  the  strain 
rate  tensor  but  not  the  magnitudes.  For  a  rigid  plastic 
solid,  the  velocities  and  strain  rates  in  problems  in  which 
inertial  forces  may  be  neglected  can  be  referred  to  a  suitable 
monotonically  varying  parameter  of  the  problem  under  considera¬ 
tion.  This  is  true  for  both  perfectly  plastic  solids  and 
plastic  workhardening  solids.  The  Levy-Mises  relation  for  a 
rigid  perfectly  plastic  solid  may  be  written  in  the  form  (13) 


k 


(2.11) 


(if  J  =  k)  , 


and  this  clearly  indicates  that  only  the  ratios  of  the  compo¬ 
nents  of  the  strain  rate  tensor  are  defined  by  the  stress 
during  plastic  deformation.  The  constitutive  relationships 
given  by  Hill  (14)  for  a  workhardening  solid  also  indicate  this. 


relation  (2.11)  to  form  the  relation  for  Bingham  solids  (2.6), 
the  magnitudes  of  the  strain  rates  and  not  just  their  ratios 
are  defined  by  the  state  of  stress  when  J  >  k.  In  a  visco¬ 
plastic  problem,  all  rates  must  be  with  respect  to  time. 
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2.4  Linearization  of  the  Constitutive  Relation 


The  linearization  of  the  constitutive  relation  of 
the  Bingham  solid  by  Prager  (5)  is  summarized  in  this  section. 
Prager  noted  that  the  constitutive  relation  in  (2.9)  and  (2.10) 
can  be  put  in  the  form, 


e .  . 


c)F 

5o.  . 


(2.12) 


where 


F  -  s/2  (J  -  k)  =  s± j s± j  "  V 2  k  (2.13) 


and 


for  F  <T  0 
for  F  ^  0 


(2.14) 


In  (2.13), the  yield  function,  F,  is  a  single  non-linear 
function.  The  linearization  involves  the  approximation  of 
this  function  by  a  set  of  m  linear  functions, 
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The  number  m  may  be  chosen  according 


to  the  degree  of  accuracy  desired.  With  the  functions  L 
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in  place  of  F  in  (2.12),  and  tu  for  the  viscosity  coefficient:, 
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Since 
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(2.16) 


then 
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q=l 
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(2.17) 


Equation  (2.17)  then  gives  the  linearized  constitutive  relation. 


The  viscosity  coefficients  r\  and  r]L  are  not,  in 
general,  the  same;  however,  they  may  be  made  the  same  by  mul¬ 


tiplying  the  functions  L 
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by  an  appropriate  constant  fac¬ 


tor.  For  example,  referring  to  Fig.  2.2,  the  Tresca  yield 
criterion  gives  the  following  linear  functions: 


Line 

AB 

"  Ld) 

=  a±  - 

°2  ' 

2k 

Line 

BC 

"  L<2) 

II 

Q 

I-1 

1 

°3  - 

2k 

Line 

CD 

"  L ( 3) 

=  a2  - 

°3  " 

2k 

Line 

DE 

-  L(4) 

=  a  2  - 

°1  " 

2k 

Line 

EF 

"  L(5) 

=  °3  - 

°1  - 

2k 

Line 

AF 

'  L(6) 

II 

Q 

OJ 

1 

°2  - 

2k 

(2.18) 


OQ  ctid 

11 


Fig.  2.2  Yield  Locus  in  Plane  a0  =  constant 
for  Tresca  Yield  Criterion 

Then,  from  (2.16),  the  stress-strain  rate  relationships  are 
as  follows:  If  the  stress  point  is  in 
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and  so  forth  for  the  other  regions. 


Now  if  k  approaches  zero,  the  following  relationships  are 
obtained; 
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From  comparing  the  linear  relationships  with  the  non-linear 
relationships  having  k— ►O,  it  follows  that 


rjL  =  3  r|  .  (2.19) 

This  difference  in  the  viscosity  coefficients  may  be  elimina¬ 
ted  by  multiplying  the  linear  functions  (2.18)  by  the  constant 
factor 

v/3 


2 . 5  The  Boundary  Value  Problem 

The  statement  of  the  boundary  value  problem  is  as 
follows:  A  visco-plastic  material  occupies  the  volume  V  and 

is  bounded  by  surface  S.  The  body  force  per  unit  volume  F 
is  given  in  V.  Surface  traction  T^  is  prescribed  on  surface 
ST,  a  portion  of  S,  and  the  velocity  v^  is  prescribed  on 
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surface  ,  another  portion  of  S .  ST  and  Sy  may  overlap, 

provided  on  the  overlapped  surface  the  prescribed  surface 

tractions  and  velocity  complement  each  other.  For  example  the 

normal  component  of  stress  and  the  tangential  component  of 

velocity  may  be  prescribed  or  vice-versa.  The  boundary  value 

problem  is  then,  for  given  and  the  prescribed  data  on  the 

surface,  to  determine  the  velocity  field  v^  and  the  stress 

field  a. .  throughout  the  volume  V.  The  stress  field  a  . 

D  D 

must  satisfy  the  equilibrium  equations 

0. .  .  +  F .  =  0  (2.20) 

ij  ,j  i 

assuming  that  inertia  effects  may  be  neglected,  and  also  it 
must  satisfy  the  boundary  condition 


a. .n.  =  T.  on  S_  .  (2.21) 

ID  D  i  T 

The  corresponding  velocity  field  v^  is  related  by  the 
stress-strain  rate  relation.  For  the  incompressible  Bingham 
solid,  the  relation  is 
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For  the  case  of  incompressible  solids,  must  also  satisfy 

the  incompressibility  relation, 


v 


1,1 


0  . 


(2.23) 


Also,  must  be  equal  to  the  prescribed  value  at  the  surface 

S^.  The  solution  of  this  boundary  value  problem  was  proved 
to  be  unique  by  Prager  (7) . 


2.6  Extremum  Principles 

The  utilization  of  variational  principles  to  obtain 
solutions  to  physical  problems  has  been  a  powerful  mathematical 
method  in  classical  and  modern  mechanics.  The  method  depends 
on  the  formulation  of  extremum  principles.  For  visco-plasticity, 
Prager  (7)  derived  two  extremum  principles  for  an  incompressible 
Bingham  solid,  under  the  boundary  condition  that  surface  trac¬ 
tions  and  velocity  are  prescribed  on  ST  and  Sv,  respectively, 
and  and  ST  are  not  zero.  Earlier,  Ilyushin  (15)  has 

derived  one  of  these  principles  assuming  =  0. 

In  order  to  facilitate  the  stating  of  the  two  extremum 
principles,  two  types  of  admissible  field  shall  first  be  defined. 
A  kinematically  admissible  velocity  field  for  a  visco-plastic 
problem  is  a  continuous  velocity  field  satisfying  the  boundary 
condition  on  and  the  incompressibility  condition.  The 
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strain  rate  field  derived  from  this  velocity  field  furnishes 
a  stress  field.  In  general,  this  stress  field  will  not  satisfy 
the  equilibrium  equations  and  the  stress  boundary  condition.  A 
statically  admissible  stress  field  is  a  stress  field  that  satis¬ 
fies  the  equilibrium  equations  and  the  stress  boundary  condition 
on  S,^ ,  The  strain  rate  field  derived  from  this  stress  field, 
in  general,  does  not  satisfy  the  compatibility  equations  and 
cannot  be  derived  from  a  velocity  field. 


The  minimum  principle  for  velocity  fields  states  that 


among  all  kinematically  admissible  velocity  fields,  the  actual 
velocity  field  minimizes  the  functional 


(2.24) 


This  principle  holds  only  for  problem  with  f  0.  Functional 

(2.24)  is  similar  to  the  potential  energy  in  classical  elasti¬ 
city  and  plasticity?  however,  the  first  integral  does  not  rep¬ 
resent  the  total  rate  of  energy  dissipation,  which  is 


r 


(2.25) 


The  maximum  principle  for  stress  fields  states  that 


among  all  statically  admissible  stress  fields,  the  actual  stress 
field  maximizes  the  functional 
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n- 


K  =  2 


J 


T  .  v .  dS 
1  l 


V 


1_ 

4T) 


r 


v 


(  |J  -  k|  +  J  -  k)2  dV  .  (2,26) 


There  is  no  closed  comparison  between  this  functional  and  the 
complementary  energy  in  classical  elasticity  and  plasticity. 
This  principle  holds  only  for  problem  with  f  0,  These 
two  extremum  principles  can  also  be  shown  to  be  true  hen 
ST  and  S^7  overlap. 
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CHAPTER  III 


EXTENSIONS  OF  THE  EXTREMUM  PRINCIPLES  FOR  A  BINGHAM  SOLID 


Upper  and  lower  bounds  for  the  rate  of  work  done  by 


the  unknown  surface  tractions  are  easily  obtained  from  the 
extremum  principles  for  a  rigid  plastic  solid.  In  visco¬ 
plasticity,  the  extremum  principles  stated  in  Section  2.6 
are  not  readily  converted  to  the  form  for  obtaining  upper  and 
lower  bounds  for  the  rate  of  work  done  by  the  unknown  surface 
tractions.  In  this  chapter,  extensions  of  the  extremum  princi¬ 
ples  are  presented  from  which  the  bounds  may  be  obtained. 

3.1  An  Extension  of  the  Minimum  Principle  for  Velocity  Fields 

Problems  for  which  the  tractions  on  ST  and  the  body 

* 

force  F^  are  zero  are  considered.  If  v^  is  a  kinemati¬ 
cally  admissible  velocity  field,  v^  is  the  solution  velocity 


* 


field  of  a  problem  of  visco-plastic  flow,  and  I  and  I  are 
the  corresponding  invariants,  then  from  (2.24) 


Since  body  force  and  prescribed  surface  tractions  are  zero  and 


it  follows 


til  mo 


,CO;,li  1 


-  .• 


i 


(  k  .  \C  ' 
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from  (2.25)  and  the  principle  of  virtual  work+  that 


J 


V 


TividS  =  I  (rjl  +  kl)dv  . 
V 


(3.2) 


Combining  (3.1)  and  (3.2)  gives 


r 


(t)I*2  +  2kI*^dV  >  |  kldV  +  j  T^dS .  (3.3) 


V 


V 


r 


V 


The  term 


V 


kldV  represents  the  rate  of  plastic  energy  dissi¬ 


pation.  Hill  (14)  and  Prager  and  Hodge  (13)  have  given  the 
extremum  principles  for  the  rigid  perfectly  plastic  solid  and 
the  extremum  principle  concerned  with  velocity  fields  states 


that  among  all  kinematically  admissible  velocity  fields 
v^°  the  actual  one  minimizes 


r 


B°  = 


kl  dV  + 


[v°J 


dS 


T  .  v .  dS 

l  l 


(3.4) 


V 


where  fv°J  is  the  magnitude  of  velocity  discontinuities  in 


the  field  v^  and  are  the  surfaces  of  discontinuity. 


+ 


Let  v^  be  a  continuous  velocity  field  and  be  an  entirely 

unrelated  stress  field  satisfying  the  equilibrium  equations, 
further  e_. ..  =  ~:(v_.  +  v..  .)  and  then  the 


=  ”  (v .  .  +  v  .  . ) 

13  2  1,3  3,1 

principle  of  virtual  work  states  that 


a .  . e .  . dV 
13  13 


T .  v .  dS 
1  1 


+ 


F  .  v .  dV 
1  1 


A  kinematically  admissible  velocity  field  for  a  rigid 
perfectly  plastic  solid  may  have  velocity  discontinuities. 


\  :  j 

. 


V 
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Inequality  (3.3)  can  be  modified  to  give  an  upper 

bound  for  T^v^dS  if  the  velocity  solution  to  the  corres- 

SV 

ponding  problem  for  a  rigid  perfectly  plastic  solid  is  known. 

* 

If  vi  is  chosen  as  the  velocity  field  of  the  rigid  perfectly 

* 

plastic  solution  and  if  v^  has  no  surfaces  of  discontinuity 
and  v^  is  the  velocity  field  of  the  visco-plastic  solution, 
it  follows  from  (3.4)  that 


dV  < 


kldV  , 


V 


V 


(3.5) 


since  the  tractions  on  ST  are  zero. 

The  following  inequality  is  obtained  by  combining  inequalities 
(3.3)  and  (3.5) . 


^  O  <^> 

(t^I  +  kl  )dV  >  /  T_.v_.dS 


i  i 


(3.6) 


V 


V 


Inequality  (3.6)  gives  an  upper  bound  for  /  T.v.dS.  The  left 

vj  S  1 

hand  side  of  (3.6)  must  of  course  be  evaluated  from  the  velocity 
field  of  a  complete  rigid  perfectly  plastic  solution  with  no 
velocity  discontinuities.  When  the  viscous  effect  is  small 
compared  with  the  yield  effect,  that  is  when  <<  1,  where 
V  is  a  characteristic  velocity  and  L  a  characteristic  length, 
it  is  reasonable  to  suppose  that  inequality  (3.6)  gives  a  close 
upper  bound. 


o  i  d  f/OA 


. 


:  •  .  J  :  . 


.  !  r3<  '*  i£  ■  >  (0  ')  on 


•.  '  '  ' 
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Many  important  rigid  perfectly  plastic  solutions 
involve  velocity  discontinuities.  If  the  velocity  field 
v^"*"  of  the  corresponding  complete  rigid  perfectly  plastic 
solution  has  velocity  discontinuities  then  it  follows  from 


where  as  before  I  is  obtained  from  the  unique  visco-plastic 
velocity  field.  The  following  inequality  is  obtained  by  com¬ 
bining  inequalities  (3.3)  and  (3.7). 


[  (t}I*2+  2kI*)dV  - 

J  kI+dV  +  \ 

["  k  [v+J  dS 

V 

__  J  v 

S, 

> 


T  .  v  .  dS  , 

l  l  * 


(3.8) 


V 


where  I  is  obtained  from  any  velocity  field  that  is  kine¬ 
matically  admissible  for  the  rigid  visco-plastic  boundary  value 

problem.  Inequality  (3.8)  then  gives  upper  bound  for  /  T . v , dS 

U  c  li 

when  the  corresponding  rigid  perfectly  plastic  solution 
has  velocity  discontinuities. 


V 


3.2  An  Extension  of  the  Maximum  Principle  for  Stresses 

This  extension  is  also  for  problems  for  which  the 

surface  tractions  on  and  body  force  are  zero.  If 

* 

a. .  is  a  statically  admissible  stress  field  ,  c. .  is  the 
lj  13 

unique  stress  field  solution  of  the  visco-plastic  problem,  and 
J  and  J  are  the  corresponding  stress  invariants,  then  from 


(2.26) 
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I  Tivids  -  jr  /  [|j-k|+  J-k]2dV>2  f  Ti*vidS 
sv  Jv  \ 


1 

4tj 


V 

f  [l  J*-k|+  J*-k]  2dV.  (3.9) 
J  V 


It  is  now  required  that  be  chosen  such  that  J  ^  k. 

Then,  it  follows  from  (3  .,9)  that 


r  ,  \  c  * 

TividS  >  /  Ti  v^dS 

S  J  s 

°V 


2t)  /  <J*-k)2dV  +  ^ 

V  '  V 


(J-k)^dV,  (3.10) 


where  is  the  volume  of  the  deforming  region  and  the 

equality  sign  is  included  only  if  the  problem  considered  is 
fully  deforming.  From  (2.8) 


k  =  T]I 


if  1^0. 


(3.11) 


Substitution  of  (3.11)  in  (3.10)  gives 


x  /  *  1 

T.v.dS>  T.  v .  dS  - 


~  ,  (J*-k)  2dV  +  -3  /  I“dV, 

2r>  Jv  27vd 


1  1 


1  1 


V 


V 


(3.12) 


The  integral  J  TjI~dV  represents  the  rate  of  viscous  energy 

vd 

dissipation  of  the  deforming  region.  From  theory  of  slow 
incompressible  Newtonian  flow,  the  extremum  principle  for  the 
velocity  field  states  that  among  all  kinematically  admissible 
velocity  fields  v^  the  velocity  solution  field  minimizes 
the  functional  (16) 


_  1  |  .+2,,. 

K  =  -^  I  r)I  dV 

2  V 


T.  v.+dS. 

l  l 


(3.13) 
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Inequality  (3.12)  can  be  modified  to  give  a  lower  bound  for 

T^v^dS  if  the  velocity  solution  to  the  corresponding  slow 
incompressible  Newtonian  flow  is  known.  Let  v^°  be  the 
velocity  solution  of  the  Newtonian  flow  problem?  and  since  the 
velocity  solution,  v^  ,  of  the  visco-plastic  problem  is 
kinematically  admissible  for  the  corresponding  Newtonian  flow 
problem  and  the  surface  traction  on  ST  are  zero,  it  follows 
from  (3.13)  that 


J  I2dV  >jjll°2dV. 


(3.14) 


where  1°  is  obtained  from  v.° 

i 


In  the  integral  on  the  left 
hand  side,  the  integration  over  the  non-deforming  region  is 
zero.  Therefore, 


-2  I  I2dV>^ 
z  j  v  ^ 

d 


r 


I°2dV  . 


(3.15) 


V 


The  following  inequality  is  obtained  by  adding  the  respective 

sides  of  (3.12)  and  (3.15). 

r  r  r  r 


T  .  v .  dS  >  I  T  .  *v .  dS  -  ”  (J*-  k)  2dV  +  --J 

-Sv  1  1  Jsv  1  1  2^v 


I°2dV  (3.16) 


V 


T  and  J  are  obtained  from  any  statically  admissible  stress 

-j- 

field  with  J  ^>k,  and  1°  is  obtained  from  the  velocity  solu¬ 
tion  of  the  corresponding  slow  Newtonian  flow  problem.  In¬ 
equality  (3.16)  then  gives  a  lower  bound  for  j  Tivi<*S* 


V 
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A  statically  admissible  stress  field  a. .  with  J  >  k 

-L  3 

may  be  obtained  from  a  kinematically  admissible  strain  rate 
field  e-Lj+>  which  is  nowhere  zero,  by  solving  the  non- 
homogeneous  perfectly  plastic  problem  with  distribution  of 
shear  yield  stress 


k'  =  k  +  n  /  2e .  ."**  e .  .**" 
'J  ID  ID 


(3.17) 


That  is 


*  /  1  *  *  / 

J  =  /  -z  s  .  .  s  .  .  -  k+,n/2e 

J  2  id  ID  ‘J 


kl  ekl 


+ 


(3.18) 


It  should  be  noted  that  the  strain  rate  field  e^j 

* 

associated  with  according  to  the  constitutive  relation¬ 

ships  (2.10)  is  not  kinematically  admissible  and  is  different 
from  however,  their  invariants  are  equal,  that  is, 

J  2e .  ."*"e .  =  2e.  .  e.  .  „  With  (3.18)  substituted, 

ID  ID  'J  !D  !D 

(3.16)  becomes 


T  .  v .  dS  > 

l  l 


V 


sv 


T  .  v  .  dS 

l  l 


H 

2 


I+2dV 


V 


+  * 


I°2dV, 


(3.19) 


One  lower  bound  is  easily  obtained  by  letting  e^j  be  equal 
to  the  strain  rate  solution  eij°  ° f  the  corresponding  slow 
Newtonian  problem.  Then  the  lower  bound  is  given  by 


:  I 
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* 

T.  v.dS 
l  l 


(3.20) 


where  is  associated  with  ck  ^  which  is  obtained  from 

solving  the  non-homogeneous  perfectly  plastic  problem  with  shear 
yield  stress  distribution  equal  to  k  +  J  2e^j°e^j°  . 


Another  possibly  better  lower  bound  may  be  obtained 
by  letting  eij+  ke  a  kinematically  admissible  strain  rate 
field  that  is  nowhere  zero  and  that  depends  on  some  unknown 
parameters  and  finding  the  parameters  that  maximizes  the  brac¬ 
keted  expression  in  (3.19).  This  maximizing  process  can  be 
accomplished  only  if  the  statically  admissible  stress  field 

■J* 

obtained  from  e. .  does  not  have  any  stress  discontinuities. 

13 


An  application  of  the  extensions  will  be  presented 


later  in  Chapter  IV. 
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CHAPTER  IV 


VISCO-PLASTIC  TORSION 


In  this  chapter,  the  problem  of  pure  torsion  of  a 


right  prismatic  bar  is  extended  to  a  homogeneous  isotropic 
visco-plastic  material,  the  Bingham  solid.  The  initial  yield¬ 
ing  of  the  solid  is  assumed  to  be  governed  by  the  Mises  yield 
criterion.  Unlike  the  plate  problem  solved  by  Appleby  and 
Prager  (10) ,  this  problem  is  not  simplified  by  assuming  the 
Tresca  yield  criterion  or  any  piece-wise  linear  yield  criterion 
instead  of  the  Mises  yield  criterion.  For  the  torsion  prob¬ 
lem,  all  isotropic  yield  criteria  that  do  not  exhibit  any 


Bauschinger  effect  reduce  to 


and  the 


yield  locus  in  the  (s  ,  s  )  plane  is  a  circle,  as  shown 

xz  yz 


in  Fig.  4.1.  Any  attempt  to  linearize  this  locus  results 


s 


yz 


s 


xz 


Fig.  4.1  Yield  Locus  for  Pure  Torsion 
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in  an  anisotropic  yield  criterion,  since  a  chord  drawn  on  the 
circle  has  a  distance  from  the  origin  smaller  than  k  and, 
consequently,  a  yield  stress  smaller  than  k  is  obtained  for 
different  directions  of  the  stress  vector. 

Approximate  solutions  are  obtained  for  the  torsion 
problem  by  a  semi-inverse  procedure,  an  approximating  procedure, 
and  the  Ritz 1 s  method.  The  torques  evaluated  from  these  approxi¬ 
mate  solutions  are  then  compared  with  the  upper  and  lower 
bound  torques  obtained  from  the  extensions  of  the  extremum 
principles  presented  in  Chapter  III. 

4.1  Statement  of  Problem 

The  z-axis  is  taken  parallel  to  a  generator  of  the 
bar  and  one  end  is  fixed  in  the  (x,y)  plane  as  shown  in 
Fig.  4.2.  The  statement  of  the  problem  is  as  follows: 


Fig.  4.2  Right  Prismatic  Bar  under  Pure  Torsion 
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The  bar  is  twisted  by  a  pure  torque  in  the  z-direction  and 
is  assumed  to  be  rigid  until  the  torque  has  exceeded  the  fully 
plastic  torque.  When  it  starts  to  deform,  the  rate  of  twist 
per  unit  length  ©  is  a  constant  and  the  ends  are  free  to  warp. 

4.2  Semi-Inverse  Procedure 

A  semi-inverse  procedure,  similar  to  Saint  Venant ' s 
procedure  for  the  elastic  torsion  problem,  is  used  to  obtain 
two  formulations  of  the  problem,  one  in  terms  of  the  warping 
function  and  the  other  in  terms  of  a  stress  function.  An 
approximate  solution  to  the  problem  is  obtained  by  applying  a 
perturbation  technique  to  the  formulation  in  terms  of  the  stress 
function. 


For  the  semi-inverse  procedure  in  terms  of  warping 
function,  the  following  velocity  components  are  assumed: 


u  =  -  0zy 
V  =  0ZX 

W  =  w  (x,y ,  0) 


> 


j 


(4.1) 


where  w  is  the  unknown  warping  function.  The  corresponding 
strain  rate  components  are 


•  ■ 
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xz 


yz 


e  =  e  =  e 
yy  zz  xy 

l,dw  A  \ 

2  -  0y) 

+  0x)  . 


=  0 


> 


(4.2) 


and 


I  = 


i 


0  ,  2 ,  2 
2  |  e  +  e 
xz  yz 


-  ey)  2+  (|s  +  ex) 2 


1 

2 


(4.3) 


The  strain  rates  e  and  e  are  related  by  the 

xz  yz 

compatibility  equation 


<$e  de 

xz  yz 

c3y  dx 


-  0  . 


(4.4) 


Substitution  of  the  strain  rates  (4.2)  in  the  constitutive 
relationships  (2.6)  gives 


(4.5) 


These  stress  components  must  satisfy  the  equilibrium  equation 


bs 


xz 


ds 


'X 


+ 


Sy 


yz  =  o 


(4.6) 
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consequently,  the  equilibrium  equation  in  terms  of  the  warping 
function  is 


2  ,  k 

T]V  W  i - ^ 

I 


2  b  2w 


>  2 
o  w 

SxSy 


+ 


2  d2w 


0  . 


(4.7) 


This  is  a  quasi-linear  elliptic  partial  differential  equation. 
Since  the  lateral  boundary  is  stress  free,  the  boundary 
condition  is 


s  n  +  s  n  =0 
xz  x  yz  y 


(4.8) 


where  nx  and  n  are  the  components  of  the  outward  unit 
normal  vector  of  the  boundary  curve.  From  (4.5) ,  this  becomes 


(4.9) 


Since  the  term  in  the  first  bracket  is  positive  definite,  the 
boundary  condition  to  be  satisfied  by  w  is 


or 


4^  n  4-  4^  n  =  0  (yn  -  xn  ) 

ox  x  oy  y  x  y 

ss  =  e  (ynx  -  my) 


(4.10) 


(4.11) 


at  the  boundary. 
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Therefore,  the  formulated  boundary  value  problem  in  terms  of 
warping  function  is  equation  (4.7)  subject  to  the  boundary 
condition  (4.11).  With  this  formulated  problem,  stress  field 
(4.5)  can  be  shown  to  be  equivalent  to  a  pure  torque. 


The  formulation  in  terms  of  the  stress  function 
involves  the  definition  of  a  stress  function  U(x,y) ,  which  is 
continuous  and  has  piece-wise  continuous  first  derivatives, 
such  that 


c)U 

sxz  “  3^ 


(4.12) 


(4.13) 


These  stress  components  automatically  satisfies  equilibrium 
equation  (4.6).  The  stress  invariant  and  corresponding  strain 
rates  are  then 


J  = 


e 

xz 


1 

2 


(4.14) 

Bu 

3y 

(4.15) 

dU 

3^  * 

(4.16) 
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These  strain  rate  components  must  satisfy  the  compatibility 
equation  (4.4).  After  substitution  of  (4.15)  and  (4.16)  in 
(4.4)  and  rearranging, 


?  v 
V  U  - 


2  d2U  „  du  du  d2U  .  ,dux  2  d2U 


-  2 


d"x  ay  cbSy  ^  2 


The  boundary  condition  to  be  satisfied  is  obtained 


thus 


-  2^9.  (4.17) 
from  (4.8) ; 


du  du 

—  n  -  ^ —  n 

ay  x  ox  y 


0  . 


(4.18) 


Since 


n 

x 


dy 

ds 


(4.19) 


and 


n 

Y 


dx 

ds 


(4.20) 


where  s  is  the  length  measured  along  the  boundary  curve, 
equation  (4.18)  becomes 


or 


du  dy  ,dudx_dU_r. 
c Yy  ds  ds  ds 

U  =  constant  on  the  boundary. 


(4.21) 

(4.22) 


For  simply-connected  regions,  the  constant  can  be  taken  as 
zero;  therefore, 

U  =  0  on  the  boundary.  (4.23) 

Then,  the  formulated  boundary  value  problem  in  terms  of  stress 
function  is  equation  (4.17)  subject  to  the  boundary  condition 
(4.23)  or  (4.22) . 
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4.3  Visco-Plastic  Torsion  of  a  Prismatic  Bar  of  Square 

Cross-Section 

The  problem  is  specialized  in  this  section  for  a  bar 
with  a  square  square  cross-section.  Referring  to  Fig.  4.3, 
the  solutions  of  the  corresponding  elastic  and  plastic  problems 
have  symmetry  about  the  x  and  y  axes,  and  in  each  quadrant, 
the  warping  function  and  warping  velocity  are  anti- symmetric 
about  the  diagonal  and  the  stress  trajectories  are  symmetric 
about  the  diagonal  (13)  (17) .  It  is  reasonable  to  assume  that 
the  solution  of  the  visco-plastic  problem  also  possesses  these 
symmetry  properties.  Therefore,  it  is  possible  to  consider 
only  one-eighth  of  the  region.  For  the  rest  of  this  chapter, 
the  triangular  region  AOB  is  being  considered  only. 


Fig.  4.3  Square  Cross-Section  of  the  Bar 
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The  boundary  conditions  derived  in  Section  4.2  can 
now  be  written  for  region  AOB .  For  the  warping  function 
formulation,  the  boundary  condition  (4.11)  becomes 

^  -  Ox  at  y  =  a  .  (4.24) 

To  ensure  that  the  warping  velocity  is  continuous  throughout 
the  region,  w  must  vanish  on  OA  and  OB  ,  because  of  the 
foregoing  symmetry  considerations.  Therefore, 

w  =  0  at  x  =  0  (4.25) 

and  w  =  0  at  x  =  y  .  (4.26) 


For  the  stress  function  formulation,  boundary  condition 
(4.23)  becomes 

U  =  0  at  y  =  a.  (4.27) 


Since  the  possibilities  of  stress  discontinuities  at  lines 
OA  and  OB  cannot  be  excluded,  no  condition  on  the  stress 
function  can  be  specified  on  these  lines.  However,  the 
corresponding  velocity  solution  w  derived  from  (4.15)  and 
(4.16)  in  the  form 


c)w 

dx 


-f 


''0 


02 


/*!¥ 


2*  02  J 


dU 


dU 


(4.28) 


(4.29) 


must  satisfy  (4.25)  and  (4.26). 
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4.4  Approximate  Solution  by  Perturbation  Method 


The  viscosity  effect  in  the  torsion  problem  can  be 
represented  by  a  dimensionless  parameter  €  =  .  This 

parameter  gives  the  ratio  of  the  characteristic  viscous  shear 
stress  and  plastic  shear  stress.  For  €  <<  1,  that  is  when  the 
viscosity  effect  is  small,  it  is  reasonable  to  assume  that  the 
visco-plastic  solution  is  slightly  perturbed  from  the 
perfectly  plastic  solution.  To  find  such  a  solution,  a  per¬ 
turbation  procedure  (18)  is  used  to  solve  the  boundary  value 
problem  for  stress  function,  U,  stated  in  the  previous  section. 
The  following  dimensionless  variables  are  introduced: 

(4.30) 


U 


x  = 


y  = 


€  = 


u 

ka 

x 

a 

X 

a 

nQa 

k 


Equation  (4.17)  expressed  in  terms  of  these  dimensionless 

* 

variables  becomes  (with  the  bars  dropped  ) 


<i>2  -  ‘I?'2 


3. 

2 


(V  U  +  2e) 


,dux  2 


(4.31) 


„  dU  dU  d2U  ,  2 

2  53E  57  S37  + 


afu 

Sy2 
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Refer  to  nomenclature  on  p.  x. 
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To  proceed  with  the  method,  the  solution  is  assumed  to  be  of 
the  series  form 

U  =  U  +  eU,  +  e2U0  +  .  .  .  .  (4.32) 

o  1  2 

Upon  substituting  in  (4.31),  using  the  binomial  expansion  for 
the  power  terms  and  collecting  terms  of  the  same  order  of  e, 
the  following  equations  are  obtained: 


For  zeroth  order  of  € 


3  0 

dU  2 

a2u 

du 

du 

32u 

du  „ 

32u  1 

vUo 
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For  first  order  of  e 


3  du  0 

d2u. 
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v  u 
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du 

du 

S2u 
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32u  1 
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dU 
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ox 

2r— 2 

dy 

dxdy 

2\ — — 
ox 

dy2 

du 

r\ 

dU 

32U 

dU 

r\ 

a2u  1 

o 

dU 

•\ —  + 

dy 

2-^ — — 
ox 

dxdy 

K) 

0> 

•C 

( 

N  2 

ox 

dy 

+  2 


VU, 


=  0 


(4.34) 


€  . 


.33) 


And  so  forth  for  higher  order  of 
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In  (4.33)  and  (4.34) , 


(4.35) 


The  boundary  conditions  for  the  above  equations  are  obtained  by 
substituting  series  (4.32)  into  the  dimensionless  form  of  (4.27). 
The  resulting  boundary  conditions  are 


uo  =  ui  =  u2  =  •  •  •  =  0  at  y  =  1.  (4.36) 

The  fully  plastic  stress  function  solution  for  the  corresponding 

perfectly  plastic  problem  is  U  =  1  -  y.  This  satisfies  the 

P 

zeroth  order  equation  (4.33).  Therefore,  it  is  reasonable  to 
take  this  as  the  zeroth  order  solution.  Thus, 

Uq  =  1  -  y.  (4.37) 

With  Uq  known,  can  then  be  found  from  (4.34).  Many 

terms  vanish  when  (4.37)  is  substituted  into  (4.34)  .  The  final 
equation  for  is 

d2U, 

■ - -  +2=0.  (4.38) 

8y2 

The  solution  for  this  equation  is 

U1  =  -  y2  +  y  f(x)  +  q(x)  (4.39) 

where  f(x)  and  q(x)  are  as  yet  unknown  functions  of  x. 
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From  boundary  condition  (4.36), 
q(x)  =  1  -  f (x) . 


(4.40) 


Therefore , 

Uj_  =  1  -  Y2  +  f  (X)  (y-1)  .  (4.41) 

The  function  f(x)  must  now  be  determined  from  the  velocity 

boundary  conditions  (4.25)  and  (4.26)  .  In  order  to  do  this, 

the  velocity  field  corresponding  to  the  stress  function 

U  =  U  +  eU.  must  be  obtained  from  the  dimensionless  form  of 
o  1 


(4.28)  and  (4.29)  written  below  with  w  = 


w 


9a‘ 


(Bars  actually 


omitted. ) 


dw  _  1 

^  -  y _  i 


i  - 


+  (th2 


'dy 


dU 

5? 


(4.42) 


dw  ,  _  1 

3 - r  X  =  -  — 

oy  e 


1  - 


2  -  $>2 


U  =  l-  y+  e^l-y2+  f  (x)  (y  -  1)J 


au 

^x 


xz 


yz 


au 


^  =  -1  +  e  [  ~2y  +  f  (x)] 


=  -  |^  =  -e  f  '  (x)  (y  -  1) 


(4.43) 

(4.44) 

(4.45) 

(4.46) 


where  s  and  s  are  in  dimensionless  form  (i.e.,  divided  by  k) 
xz  yz  ’ 
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M 


2  ,  ,du*  2 

y 


+  w 


■ 

r  i  ?  r  i 

=  J 

1  +  2e  2y  -  f(x)  +  €  2y-f(x) 

2+  € 2 

f ' (x) (y-1) 

L  J  L  J 

L-  J 

nj 


1  +  6  [2y  -  f(x)] 


(4.47) 


Therefore,  upon  substitution  of  (4.44)  to  (4.47)  and  using 
binomial  expansion,  (4.42)  and  (4.43)  become 


dw 


^  -  y  =  -  2y  +  f (x) 


dw  ,  _ 

3?  +  x  "  e 


f '  (x)  (y  -  1)  \_2y  -  f  (x)_ 


From  (4.48) , 


w  =  -xy  4-  1  f(x)dx  +  h(y)  , 


(4.48) 


(4.49) 


(4.50) 


where  h(y)  is  an  as  yet  unknown  function  of  y.  From 
boundary  condition  (4.25),  that  is,  w  =  0  at  x  =  0, 


h(y)  -  constant  =- 


f (x) dx 


J  x=0 


Therefore, 


r 


X 


w  =  -xy  + 


J 


f  ( z )  dz  . 


(4.51) 


From  boundary  condition  (4.26) ,  that  is  w  =  0  at  x  =  y, 


x 


f  (z)  dz  =  x‘ 


and 


f(x)  =  2x„ 


(4.52) 
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(4.52)  will  not  satisfy  equation  (4.49)  exactly,  but  it  is  satis¬ 
fied  to  the  order  of  approximation  of  the  perturbation  method, 
since  the  left  hand  side  is  zero  and  the  right  hand  side  has  a 
maximum  error  equal  to  €  and  a  root  mean  square  error  of 
0.298e.  Therefore,  for  €<<1,  the  approximate  solution  of  the 
problem  may  be  taken  as 


s 


xz 


1  -  2e  (y  -  x) 


> 


(4.53) 


s 


yz 


=  2e  (1  -  y)  . 


The  stress  trajectories  or  constant  U  lines  are  shown  in 
Fig.  4.4.  The  torque  in  dimensionless  form  is 


(syzx  "  sxzy)  dxdY  • 


(4.54) 


With  the  approximate  solution  (4.53) 


T  =  -|  (1  +  e)  . 


(4.55) 
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Fig.  4.4  Approximate  Stress  Trajectories  from 

Perturbation  Method 
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4.5  Approximate  Solutions  -  Method  (a) 

In  order  to  have  other  approximate  solutions  to  com¬ 
pare  with  the  solution  obtained  by  the  perturbation  method,  two 
other  approximate  methods  were  applied.  The  principle  of  the 
first  approximate  method,  henceforth  referred  to  as  Method  (a) , 

is  as  follows:  From  (4.5)  or  (2.6),  the  magnitude  of  the  shear 

* 

stress  at  each  point  of  the  cross-section  is 


(4.56) 


From  this  equation,  it  can  be  deduced  that  the  visco-plastic 
torsion  problem  has  the  same  stress  solution  as  the  corresponding 
problem  of  perfectly  plastic  torsion  of  a  bar  with  non-uniform 
hardness  such  that  the  yield  stress  in  shear  has  a  distribution 


(4.57) 


k 1  =  k  +  t)I 


in  which  I  is  a  function  of  x  and  y,  and  is  derived  from 
the  velocity  solution  of  the  visco-plastic  problem.  This  analogy 
does  not  help  in  getting  an  exact  solution;  however,  an  approxi¬ 
mate  solution  may  be  obtained  by  deriving  I  from  the  corres¬ 
ponding  perfectly  plastic  velocity  solution  instead  of  from  the 

■jlf 

Equations  (4.56)  to  (4.60) ,  inclusive,  are  in  dimensional  form. 
Refer  to  nomenclature  on  p.  x  for  the  others. 

The  velocity  field  of  the  non-uniform  perfectly  plastic  problem 
is  the  same  as  that  for  the  visco-plastic  problem  to  within  an 
arbitrary  multiplier  only. 
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visco-plastic  velocity  solution  and  solving  the  inhomogeneous 
perfectly  plastic  torsion  problem,  for  which  a  method  of  solution 
is  known  (19) ,  with  yield  stress  distribution 


where 


k '  =  k  +  Tjlp  , 


c>w  .  2 

<3?  "  0y)  + 


(4.58) 


(4.59) 


and  Wp  is  the  velocity  solution  of  the  corresponding  perfectly 
plastic  problem.  For  a  square  section,  referring  to  region 
AOB  in  Fig.  4.3, 


w  =  -  6x(y  -  x) 

Jr' 


(4.60) 


In  dimensionless  form 


w  =  -x (y  -  x) 

Jl 


(4.61) 


Therefore,  yield  stress  distribution  (4.58)  in  dimensionless 
form  is 


P  =  *  1  +  el 

K  k  p 


1  +  € 


dw  0  8w  0 

<3?  -  y)  +  x) 


n  1 


=  1  +  2e  (y  -  x) 


(4.62) 


The  dimensionless  shear  stress  may  be  written  as 

s  =  -  6  cos  <t> 
xz 

s  =  p  sin  <t>  (4.63) 

yz  * 


where  $  is  the  angle  between  the  shear  stress  vector  and 
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the  direction  of  the  negative  x-axis.  Substituting  these  stress 
components  into  equilibrium  equation  (4.6)  gives 


A  ,  A 

sin  <p  ^  +  cos  <p 


=  i  (cos  $  || 


sin  0  ||) 


or 


=  -  i  (cos  $  ||  - 


sin  *  §|)  , 


(4.64) 


(4.65) 


=  -  sin  $  -  cos  <t>  ^  and  n  is  in  the 


where  ^ —  —  —  sm  ^  ^ — 
on  ox 

direction  normal  and  90°  counterclockwise  from  the  direction  of  the 


stress  trajectories,  and  is  in  the  dimensionless  form,  that  is, 
divided  by  a.  Since  the  boundary  of  the  cross-section  is  a 
stress  trajectory,  <t>  is  known  on  it.  Therefore,  (4.65)  may 
be  integrated  inward  from  the  boundary  and  normal  to  it  to  deter¬ 
mine  ^  for  the  whole  cross-sectional  region.  The  integration 
of  (4.65)  or  (4.64)  is  done  by  a  numerical  procedure. 

Equation  (4.64)  is  used  for  the  numerical  procedure, 
because  it  is  easier  to  work  with  square  grid  points  than  irregu¬ 
lar  grid  points  in  which  utilization  of  (4.65)  will  result. 

^  cj)  ,  6 

With  forward  difference  for  ^  and  backward  difference  for  , 

for  a  typical  grid  point  (i,j),  (4.64)  in  finite  difference  form  is 

tan  ^  (i,  j)  ~  +  j)  ~  Vi,  j-1) 


—  h 


^fx^  (i,  j)  ”  tan  $(i,j)  (3y^  ( i ,  j ) 


,  (4.66) 


.  ic  :  ■  :  *X  b 

-  -  XU  a  ■  ■  ) 


■ 


.  'O  -  I 


.  :  -  ■ 


.  ■  .  Of?  hi 

■■■  .  .  ’  '  :  '  -  II- 

"  •  -  V  -•  - 

' 

^urii 

■  ■  .  .. 

•  ,  . '  • ;  j  •  ,  '  .  •  -  -ryZ 

■ 

.  s.  i  ■ 


■ 


45 


where  h  is  the  grid  width.  This  equation  relates  the  values  of 
0  at  three  grid  points  (i+1,  j) ,  (i,j),  and  (i,j-l).  Referring 
to  Fig.  4.5,  the  numerical  calculation  is  started  at  the  upper 
left  hand  corner  of  region  AOB.  Since  is  known  at  grid 

points  A  and  D,  0  at  grid  point  C  can  be  determined  from 
(4.62)  and  (4.66).  By  proceeding  to  the  right,  <t>  for  all  grid 
points  on  line  CE  can  be  determined.  In  the  same  manner, 

<t>  for  the  next  grid  line  below  can  be  determined,  and  so  forth 
for  the  whole  region. 


Fig.  4.5  Typical  Grided  Region  for  Numerical  Calculation 
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Fig.  4.6 


Stress  Trajectories  from  Approximate 

Method  (a) 
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The  above  numerical  procedure  is  performed  by  Computer 

* 

Program  A  ,  which  also  calculates  the  stress  components  from 
(4.63)  and  the  dimensionless  torque  from  (4.54).  For  e  =  0.1 
and  grid  size  of  0.01  of  the  width  of  the  square,  values  of  4> 
and  stresses  were  calculated  and  tabulated  in  Table  A.l  of  the 
Appendix,  and  the  dimensionless  torque  is  2.9281.  The 
approximate  stress  trajectories  are  shown  in  Fig.  4.6. 


4.6  Application  of  Ritz’s  Method  -  Approximate  Method  (b) 


Ritz ' s  method  is  used  along  with  the  extremum  princi¬ 
ple  for  velocity  fields  to  obtain  approximate  solutions  to  the 
torsion  problem.  To  apply  the  velocity  field  extremum  principle 
to  the  torsion  problem,  it  is  assumed  that  the  rate  of  twist 
per  unit  length  6  is  prescribed  by  specifying  the  velocities 
(referring  to  Fig.  4.2)  in  dimensional  forms 


u  =  0 
v  «  0 


at  z  =  0 


u  =  -  Q'ty 


v  =  Qtx 


> 


at  z  =1  . 


The  lateral  surface  is  stress  free  and,  since  the  ends  are  free 
to  warp,  the  normal  component  of  stress  on  the  ends  is  prescribed 


* 


All  the  computer  programs  are  given  in  the  Appendix. 
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to  be  zero.  Also  the  body  forces  due  to  the  weight  of  the  material 
are  assumed  to  be  zero.  From  (2.24)  the  functional  to  be 
minimized  is  (in  dimensional  form) 


* 

H 


ie  O  *k 

(tj!  +  2k I  )  dV 


(4.67) 


where  I  is  obtained  from  a  kinematically  admissible  strain 
rate  field.  For  this  pure  torsion  problem,  dV  =  £dS ,  where 
l  is  the  length  of  the  bar.  Since  it  is  simpler  to  work  with 
dimensionless  quantities,  the  following  dimensionless  variables 

are  introduced: 

* 


— * 


H 


H 

*  4  > 

0a  k 


I 


I_ 

o 

0a 


9 


(4.68) 

_  n  _  dS 

l  =  -  ,  ds  =  2 

a  *  a 


Therefore,  the  dimensionless  form  of  (4.67)  is  (with  the 
bars  dropped) 


* 

H 


=  l 


■k  O  * 

(■€  I  +  21  )  dS 


s 


(4.69) 


where 


* 

I  = 


rdw 

kcTx~ 


-  y)  + 


.bw 


+  x) 


1 

2 


(4.70) 


Two  kinematically  admissible  velocity  fields  are  chosen. 

* 


w 

* 

w 


=  yvp  + 


(1  -  7)wv 


They  are 

(4.71) 

(4.72) 
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where  is  the  fully  plastic  velocity  solution  for  the 

corresponding  perfectly  plastic  problem;  is  the  velocity 

solution  for  the  corresponding  Newtonian  flow  problem;  and 

* 

y  is  an  unknown  multiplication  factor  to  be  chosen  so  that  H 
is  minimized.  Functions  depending  on  only  one  parameter  are 
chosen,  because  the  highly  non-linear  nature  of  the  functional 
(4.69)  would  make  the  minimizing  process  quite  difficult  for 
more  than  one  parameter.  Henceforth,  the  application  of  Ritz ' s 
method  using  velocity  field  (4.71)  will  be  referred  to  as 
Method  (b-1)  and  that  using  velocity  field  (4.72)  as  Method  (b-2) . 

To  apply  Method  (b-1) ,  velocity  field  w  for 

P 

region  AOB  from  (4.61)  is  substituted  in  (4.71)  giving 

w  =  -yx(y  -  x)  .  (4.73) 


The  functional  (4.67)  for  this  velocity  field  and  for  a  square 
section  is 


/*■ 


_  _  H  _ 

G  "  sT  -  e  j 


AOB 


( 1+7)  y  -  2yx 


+  (l-y)2x2> 


dS 


r 


+  2 


j 


AOB 


r 


n 


( 1+7)  y  -  27 x 


+  ( 1-7) 2x2 


1 

2 


dS.  (4.74) 


It  is  easier  to  find  the  value  of  7  that  minimizes  G  by 
numerically  integrating  (4.74).  For  €  =  0.1,  the  minimizing 
process  was  done  by  a  computer  program  and  found  that  G  is 
minimum  at  7  =  1.0. 
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The  dimensionless  stresses  associated  with  the  velocity  field 
(4.73)  with  y  =  1.0  are 


2e(y  -x) 


(4.75) 


and  are  equivalent  to  a  pure  torque.  With  this  stress  field, 
the  dimensionless  torque  obtained  from  (4.54)  is 


T  =  |  (1  +  |)  . 


(4.76) 


To  apply  Method  (b-2) ,  the  pure  viscous  velocity 
solution  for  slow  Newtonian  flow  in  the  torsion  problem  is 
obtained  from  replacing  twist  by  rate  of  twist  in  the  elastic 
warping  function.  Hence,  in  dimensionless  form 


32  'X-'  (-l)n 

it3  £?0  ( 2n+l)  3 


sinh  k  y 
cosh  k 

n 


sin  k  x  , 
n  * 


(4.77) 


where  k  =  ~  (2n+l) .  Then  (4.72)  becomes 
n  b 


* 

w 


-yx(y-x) 


+  (1-7) 


32 

3 


oo 


n=0 


(-l)n 

(2n+l)  3 


sinh  k  y 
cosh  k 

n 


sin  k  x 
n 


(4.78) 


The  velocity  field  (4.78)  was  substituted  in  (4.69)  and  (4.70), 
and  the  integral  evaluated  numerically  by  a  computer  program  for 
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values  0.90 -$r  7  ^  1.00  and  e  =  0.1.  Since  the  series  in 
(4.78)  converges  quite  rapidly,  it  was  approximated  only  by  the 

* 

first  eight  terms  in  the  evaluation.  The  results  showed  that  H 
is  minimum  at  y  =  0.96  within  the  accuracy  of  the  numerical 
integration.  With  y  =  0.96  substituted  into  (4.78), 


oo 


”  y  =  -0 . 96  (y  -  x)  - 


0.64 


tt  ^  (2n+l) 


/  . v  n  sinh  k  y 

■L7-J  0  cos  k  x  (4.79) 

2  cosh  k._  n 


n 


dw  ,  _  ^  0.64 

^ —  +  x  —  0 . 08x  —  - ^ — 


oo 


Jbil 


n  cosh  k  y 


n- 


t r 


/  .  /  „  ,  x  2  cosh  k 

^=o  (2n+1) 


sin  k  x 

n 


(4.80) 


The  corresponding  dimensionless  stresses  are  then 

sxz  =  (e  +  I*  (t^  ~  y)  (4-81) 

syz  =  (€  +i)(g+  x)  (4.82) 

The  dimensionless  torque  can  then  be  evaluated  by  (4.54).  The 
numerical  values  of  the  stresses  and  the  torque  may  be  obtained 
numerically  by  Computer  Program  B.  For  e  =  0.1  and  grid  size  of 
0.01  of  the  width  of  the  square,  values  of  stresses  were  evaluated 
and  tabulated  in  Table  A. 2,  and  the  dimensionless  torque  is 
equal  to  2.9293.  Approximate  stress  trajectories  are  drawn  from 
the  values  of  the  stresses  and  shown  in  Fig.  4.7. 
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Fig.  4.7 


Stress  Trajectories  from  Approximate 

Method  (b-2) 
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4.7  Upper  and  Lower  Bound  Torques 


The  extensions  of  the  extremum  principles  presented 


in  Chapter  III  are  applied  in  this  section  to  the  visco-plastic 
torsion  of  a  bar  with  square  cross-section  to  obtain  upper  and 
lower  bound  torques.  The  results  are  used  for  comparison  with 
the  other  approximate  torques  previously  obtained. 


To  apply  the  extension  of  the  minimum  principle  for 


velocity  fields  inequality  (3.6)  is  written  for  the  torsion 
problem  in  the  dimensionless  form  as  follows: 


r 


(el*2  +  I*)dS  >  j(syzx  -  sxzY)ds 


(4.83) 


where  the  integrations  are  over  the  cross-sectional  area  of  the 

bar.  The  right  hand  side  is  the  dimensionless  torque  for  the 

torsion  problem  shown  in  (4.54);  therefore,  the  left  hand  side 

* 

is  the  upper  bound  torque.  I  must  be  evaluated  from  the  per¬ 
fectly  plastic  velocity  solution,  w  .  From  (4.61), 

P 

w  =  -x(y-x)  in  region  AOB;  consequently, 

P 


I 


* 


+  <^+  x> 


dw  0  dw 


1  - 

2  2  =  2 (y-x)  .  (4.84) 


The  upper  bound  torque  is  then 

r 


★  o  * 

T  =  8  ( €  I  +  I  )  dS 

u 

J  AOB 


^  0  ^  x 


dydx  =  —■  (1+e)  .  (4.85) 
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This  is  a  close  upper  bound  torque  for  1.  For  larger 

values  of  e,  that  is  if  the  viscous  effects  are  as  much  as  or 

much  greater  than  the  plastic  yield  effects,  a  better  upper 

bound  torque  is  obtained  from  inequality  (3.8)  by  evaluating 
* 

I  from  the  velocity  solution  of  the  corresponding  Newtonian 
flow  problem  and  I+  from  that  of  the  perfectly  plastic 
problem.  For  the  torsion  problem,  inequality  (3.8)  in  the 
dimensionless  form  is 

s  x  -  s  y) dS .  (4.86) 

yz  xzJ 

in  which  the  left-hand  side  is  the  upper  bound  torque.  From 

-f-  ic 

(4.84),  I  -  2(y-x),  and  I  is  evaluated  from  the  velocity 
field  w^  in  (4.77).  With  the  first  integral  numerically 
integrated,  the  upper  bound  torque  for  large  value  of  e  is 


*  ? 

(el  ^  + 


r 


21*)  dS 


I+dS  > 


T  =  2.249  €  +  2.918.  (4.87) 

u 

The  torque  evaluated  from  the  corresponding  perfectly 
plastic  solution  is  a  trivial  lower  bound  for  the  visco-plastic 
problem.  For  the  square  section,  the  dimensionless  trivial 

f  8 

lower  bound  torque  is  8  ydydx  =  -=r  . 

J  AOB 

A  better  lower  bound  torque  may  be  obtained  from  the 
extension  of  the  maximum  principle  for  stress.  The  inequality 
(3.20)  written  for  the  torsion  problem  is  (in  dimensionless  form) 

/syzX  “  sxzy)dS  >  j(syz*X  "  SxZ*y)dS>  (4'88) 
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where  the  integrations  are  over  the  cross-sectional  area  of  the 


bar.  The  left  hand  side  is  again  the  dimensionless  torque  for 

the  torsion  problem.  Thus,  the  right  hand  side  is  a  lower  bound 

*  * 

s  and  s  are  obtained  by  solving  the  inhomogeneous 

perfectly  plastic  problem  with  shear  yield  stress  distribution 
(in  dimensionless  form) 


I 


o 


1  +  el 

dw 


<33r  -  + 


1 

2 


(4.89) 

(4.90) 


where  w^  is  the  velocity  solution  of  the  corresponding  slow 

Newtonian  flow  problem  given  in  (4.77).  The  method  of  solution 

for  the  inhomogeneous  problem  was  presented  in  Section  4.5. 

Following  the  same  numerical  procedure.  Computer  Program  C 

*  * 

calculates  the  stress  components  s  and  s  and 

*  xz  yz 

evaluates  the  integral  on  the  right  hand  side  of  (4.88) .  For 
e  =  0.10  and  grid  size  of  0.01  of  width  of  the  square,  the  dimen 
sionless  lower  bound  torque  is  2.8950. 


To  give  some  indications  of  how  close  the  upper 
and  lower  bound  torques  enclose  the  actual  torque,  upper  and 
lower  bound  torques  were  calculated  from  (4.85)  and  Computer 
Program  C  for  the  range  0.05  4  0.60,  and  from  (4.87)  and  Com¬ 

puter  Program  C  for  the  range  0.64  €412.0.  The  results  are 
tabulated  in  Tables  4.1  and  4.2. 
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Table  4.1 

Upper  and  Lower  Bound  Torques  for  €  <<  1 


e 

Upper  Bound 
Torque 

Lower  Bound 
Torque 

Difference 

Average  of  Bounds 

0.05 

2.8000 

2.7809 

0.68% 

0.10 

2.9333 

2.8950 

1.31% 

0.15 

3.0667 

3.0091 

1.91% 

0.20 

3.2000 

3.1231 

2.43% 

0.25 

3.3333 

3.2370 

2.93% 

0.30 

3.4667 

3.3509 

3.40% 

0.35 

3.6000 

3.4647 

3.83% 

0.40 

3.7333 

3.5784 

4 . 24% 

0.45 

3.8667 

3.6921 

4.62% 

0.50 

4.0000 

3.8058 

4.98% 

0.55 

4.1333 

3.9193 

5.3  2% 

0.60 

4.2667 

4.0329 

5.63% 
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Table  4.2 


Upper  and  Lower 

Bound  Torques 

for  Larqe  e 

G 

Upper  Bound 
Torque 

Lower  Bound 
Torque 

Difference 

Average  of  Bounds 

0.50 

4.0425 

3.8058 

6.03% 

1.00 

5.1670 

4.9401 

4.49% 

1.50 

6.2915 

6.0717 

3.56% 

2.00 

7.4160 

7.2015 

2.93% 

2.50 

8.5405 

8.3303 

2.49% 

3.00 

9.6650 

9.4581 

2.16% 

3.50 

10.7895 

10.5853 

1.91% 

4.00 

11.9140 

11.7120 

1.71% 

4.50 

13.0385 

12.8383 

1,55% 

5.00 

14.1630 

13.9642 

1.41% 

5.50 

15.2875 

15.0899 

1.30% 

6.00 

16.4120 

16.2153 

1.21% 

6.50 

17.5365 

17.3406 

1.12% 

7.00 

18.6610 

18.4657 

1.05% 

7.50 

19.7855 

19.5906 

0.99% 

8.00 

20.9100 

20.7154 

0.94% 

8.50 

22.0345 

21.8401 

0.89% 

9.00 

23.1590 

22.9647 

0.84% 

9.50 

24.2835 

24.0892 

0.80% 

10.00 

25.4080 

25.2136 

0.77% 

10.50 

26.5325 

26.3380 

0.74% 

11.00 

27.6570 

27.4623 

0.71% 

11.50 

28.7815 

28.5865 

0.68% 

12.00 

29.9060 

29.7107 

0.66% 
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4.8  Summary  of  Results  and  Discussion 

The  approximate  solutions  are  presented  in  terms  of 
the  stress  fields  (sx2  and  Sy2) ,  the  approximate  stress  tra¬ 
jectories,  and  the  torques  evaluated  from  the  stress  fields. 

They  are  given  in  terms  of  dimensionless  variables  so  that  the 
solutions,  especially  the  numerical  solutions,  can  be  applicable 
to  different  characteristic  dimensions,  yield  stresses,  and 
viscosity  coefficients.  The  stress  components  are  either  in  the 
analytical  form  presented  under  each  section  or  in  the  numerical 
form  presented  in  the  Appendix.  The  stress  trajectories  are 
shown  in  figures  under  each  section.  The  various  torques 
evaluated  are  summarized  in  Table  4.3.  To  recapitulate,  Method  (a) 
is  by  solving  the  corresponding  inhomogeneous  perfectly  plastic 
problem,  Method  (b-1)  is  Ritz ' s  method  using  velocity  fields 

w  =  yw  ,  and  Method  (b-2)  is  Ritz's  method  using  velocity 
P 

* 

fields  w  =Yw  +  (1  -  y)  w  . 

p  "  v 


Table  4.3 

Summary  of  Torques 
(For  e  =  0.1) 


Method  of  Solution 

Dimensionless  Torque 

T/a*k 

Upper  Bound 

2.9333 

Perturbation  Method 

2.9333 

Method  (b-2) 

2.9293 

Method  (a) 

2.9281 

Lower  Bound 

2.8950 

Method  (b-1) 

2.8000 

Trivial  Lower  Bound 

2.6667 
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One  desirable  condition  for  the  approximate  solution 
to  be  a  good  approximation  is  that  the  torque  evaluated  from  it 
lies  within  the  upper  and  lower  bound  torques  obtained  from 
the  extremum  principles.  The  solution  obtained  by  Method  (b-1) 
gives  a  torque  that  is  lower  than  the  lower  bound  torque,  and 
its  stress  field  has  s^z  =  0,  which  would  not  be  expected  in 
the  exact  solution.  Therefore,  the  solution  obtained  from 
Method  (b-1)  is  probably  not  a  good  approximate  solution. 


The  solution  obtained  by  the  perturbation  method 
gives  a  torque  that  is  equal  to  the  upper  bound  torque;  also  it 
has  SyZ  independent  of  x,  which  in  the  exact  solution  would  be 
dependent  on  x.  Therefore,  the  perturbation  solution  is  probably 
not  as  good  an  approximate  solution  as  those  obtained  by 
Method  (a)  and  Method  (b-2) . 


In  the  other  two  solutions,  s  at  points  near  the 

>  yz 

y-axis  are  smaller  than  those  at  points  near  the  diagonal. 

This  is  reasonable,  because  the  y-components  of  the  strain  rate 
due  to  the  twisting  are  larger  near  the  diagonal  than  those 
near  the  y-axis. 


In  the  solution  obtained  by  Method  (a) ,  at  points 
near  the  diagonal,  s^z  increases  monotonically  from  the  corner 
to  the  origin.  In  the  exact  solution,  s^z  near  the  diagonal 
would  probably  increase  from  zero  at  the  origin,  reach  a  maximum 
near  the  mid-point  between  the  origin  and  the  corner,  and  decrease 
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towards  the  corner.  This  is  so,  because  the  rate  effects  due  to 

the  twisting  would  be  greater  away  from  the  origin;  but  hear  the 

corner,  because  of  the  influence  of  the  geometry  of  the  section, 

s  „  would  become  smaller  again.  The  distribution  of  s  near 
yz  r  yz 

the  diagonal  in  the  solution  obtained  by  Method  (b-2)  is  similar 
to  that  one  would  expect  for  the  exact  solution.  Therefore,  for 
the  torsion  problem,  Method  (b-2)  is  probably  the  best  among  the 
four  approximate  methods  considered.  However,  this  does  not 
mean  that  this  method  would  give  good  approximate  solution  for 
other  two  dimensional  visco-plastic  problems. 

In  conclusion,  an  approximate  solution  may  be  obtained 
for  the  torsion  problem  by  the  Ritz's  method  using  a  kinematically 
admissible  velocity  field  formed  by  the  linear  combination  of 
the  velocity  solutions  of  the  corresponding  perfectly  plastic 
and  slow  Newtonian  problems,  or  by  the  approximate  method  (a) 
suggested  in  this  chapter.  A  better  approximate  solution  may 
probably  be  obtained  by  the  combination  of  these  two  methods. 
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CHAPTER  V 

EXPERIMENTAL  INVESTIGATIONS  OF  VISCO-PLASTIC  TORSION 


The  experimental  investigation  of  the  pure  torsion  of 
circular  cylinders  of  pure  lead,  which  is  known  to  exhibit  rate 
effects  at  room  temperature,  is  presented  in  this  chapter. 

The  dependence  of  the  yield  stress  of  pure  lead  on  strain  rates 
for  various  temperatures  and  strains  has  been  investigated 
by  Loizou  and  Sims  (20)  using  the  compression  test.  Some  effec¬ 
tive  stress-effective  strain  curves  obtained  by  compressing 
commercially  pure  lead  at  various  strain  rates  were  also  in¬ 
cluded  in  the  investigation  of  Kobayashi,  et  al.  (21) .  Davis 
and  Hunter  (22)  and  Johnson  (23)  have  attempted  to  assess  the 
strain  rate  sensitivity  of  materials  by  using  dynamic  indenta¬ 
tion  tests,  but  lead  was  not  considered.  All  the  investigations 
mentioned  above  did  not  consider  a  visco-plastic  model  and  the 
strain  rate  sensitivity  data  presented  are  characteristic  of 
the  experiments  performed;  consequently,  the  application  of  the 
data  to  other  problems  is  inconvenient  and  possibly  difficult. 

In  the  compression  test,  the  definition  of  stress  and  strain 
rate  is  difficult,  and  in  both  methods  of  testing  frictional 
effects  are  difficult  to  eliminate. 

The  purpose  of  this  experimental  work  is  to  investi¬ 
gate  the  characteristics  of  pure  lead  circular  cylindrical  spe¬ 
cimens  under  pure  torsion,  and  to  consider  the  possibility 
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of  assessing  the  strain  rate  sensibility  of  commercially  pure 
lead  by  representing  the  lead  by  the  Bingham  model  and  finding 
the  viscosity  coefficient.  The  advantages  of  using  the  torsion 
test  are  the  absence  of  frictional  problems,  the  ease  of  defini¬ 
tion  of  stress  and  strain,  and  the  ease  of  solution  of  the 
theoretical  problem  for  a  Bingham  solid.  The  disadvantage  is 
that  the  stresses  and  strains  are  non-uniform. 

This  experimental  work  also  investigated  briefly 
the  shortening  of  a  pure  lead  specimen  during  torsion. 

Swift  (24)  and  Hill  (25)  pointed  out  that  a  lead  cylinder 
shortens  during  torsion,  and  Swift  presented  some  data  on  the 
lengthenings  of  circular  cylinders  of  materials  like  steel, 
copper,  and  aluminum,  but  tested  only  one  lead  specimen. 

5.1  Theoretical  Consideration. 

The  solution  of  the  problem  of  pure  torsion  of  a 

* 

circular  cylindrical  bar  of  an  elastic-Bingham  solid  under 
constant  rate  of  twisting  is  easily  obtained.  The  shear  stress 
distribution  is  as  follows: 


That  is,  a  solid  that  behaves  as  a  Hookean  solid  when  the 
stress  state  is  under  the  yield  stress  and  as  a  Bingham  solid 
when  the  stress  state  exceeds  the  yield  stress. 
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t  =  |jL0r  if  0  ^  r  ( 

t  =  rjGr  +  k  if  R,  <  r  $  R 


(5.1) 

(5.2) 


where  t  is  the  magnitude  of  the  shear  stress  vector,  R  is  the 
radius  of  the  bar,  R^  is  the  radius  of  the  elastic-plastic 
boundary,  and  the  shear  stress  vector  is  directed  tangential 
to  the  circle  of  radius  r.  R^  is  the  radius  at  which  t  is 
equal  to  k;  therefore,  from  (5.1) 


(5.3) 


It  can  be  noted  from  (5.1)  and  (5.2)  that  there  is  a  stress 
discontinuity  at  R^.  The  torque  evaluated  from  (5.1)  and 
(5.2)  is 


r 

T  =  t  rds 


2. 

Tr  dr 


From  (5.3) ,  it  follows  that 


(5.4) 


T  =  7T0T) 


«4- 


+  |rrk 


R3- 

R  4^9' 


(5.5) 


For  pure  lead,  —  is  approximately  l/390;  therefore,  for  large 
angle  of  twist  such  that  (5.5)  becomes 


1*4  9  3 

T  =  —  tt0t]R  +  -^irkR 


(5.6) 
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From  this  relation,  t)  may  be  determined  by  twisting 
specimens  of  the  same  diameter  at  various  constant  rates  of 
twist  and  measuring  the  torque  after  a  certain  angle  of  twist 
such  that  «  R.  If  pure  lead  behaves  as  a  Bingham  solid, 

the  data  should  supply  a  linear  relation 

T  =  AG  +  B  . 

Therefore,  from  the  constant  A,  the  coefficient  of 
viscosity  t)  may  be  determined. 

5.2  Experimental  Apparatus 

The  experiments  consisted  essentially  of  twisting 
circular  lead  bars  at  various  rates  of  twist  and  measuring  the 
torques  during  twist.  A  variable  speed  torsion  testing  machine 
was  designed  to  perform  the  experiments.  The  machine  consisted 
of  five  main  parts: 

(1)  The  Frame  and  Twisting  Heads  -  The  main  frame 
of  the  machine  consisted  of  the  base  frame  supporting  a  base 
plate  and  four  posts  mounted  vertically  on  the  base  plate 
supporting  the  top  plate.  The  twisting  heads  were  five-inch 
three- jaw  chucks  and  were  arranged  and  mounted  as  shown  in 
Fig.  5.1.  The  shaft  of  the  moving  head  was  in  two  parts  split¬ 
ting  at  the  point  indicated  in  Fig.  5.1.  This  permitted  the 
specimen  to  shorten  or  lengthen  while  being  twisted.  A  photo¬ 
graph  of  the  frame  and  twisting  heads  is  shown  in  Fig.  5.3. 
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Worm  Gear 


Worm 


Spur  Gear 
Pinion 


Fig.  5.1  Section  View  of  the  Main  Part  of  the  Variable 

Speed  Torsion  Testing  Machine 
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Fig.  5.2  Overall  View  of  the  Torsion  Testing  Machine 


Fig.  5.3  Frame  and  Twisting  Heads 
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(2)  Drive  and  Speed  Adjustment  -  The  drive 
consisted  of  two  sets  of  gears:  a  set  of  worm  (one  tooth)  and 
worm  gear  (60  or  30  teeth) ,  and  a  set  of  spur  gear  (60  teeth) 
and  pinion  (20  teeth) .  The  gears  and  the  motor,  arranged  as 
shown  in  Fig.  5.1,  gave  gear  reductions  of  1:180  or  1:90 
depending  on  which  worm  gear  was  used. 

A  3/4  horsepower  d.c.  motor  was  used.  A  much 
larger  motor  than  required  was  chosen  in  order  to  have  a  con¬ 
stant  speed  during  testing  without  the  necessity  of  a  feedback 
control  system.  Without  feedback  control,  the  speed  was  kept 
constant  to  within  ^  1.5%  by  two  methods  of  speed  adjustment. 

One  method  was  for  high  speed  operation,  at  1200  rpm  to  1800  rpm; 
a  rheostat  was  connected  in  series  with  the  armature  and  the 
speed  was  adjusted  by  varying  the  rheostat  resistance.  The  other 
method  was  for  low  speed  operation,  at  180  rpm  to  1400  rpm; 
the  armature  was  connected  to  an  independently  adjustable  voltage 
supplied  by  a  "Variac"  variable  transformer  through  a  rectifier, 
and  the  speed  was  adjusted  by  varying  the  "Variac"  voltage. 

The  wiring,  as  shown  in  Fig.  5.6,  was  such  that  changing  from 
one  method  of  speed  adjustment  to  the  other  only  involved 
throwing  a  switch.  The  speed  adjustment  equipment  are  shown 
in  the  center  background  of  Fig.  5.4. 
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(3)  Speed  Measurement  -  The  speed  of  the  moving 
head  was  monitored  for  constancy  by  a  d.c.  generator  type  tacho¬ 
meter  connected  to  the  motor  shaft.  The  tachometer  output  was 
fed  into  the  y-channel  of  an  x-y  recorder,  which  had  the 
x-channel  in  time  sweep.  The  actual  speed  of  the  moving  head 
was  measured  by  recording  the  time  lapse  for  each  revolution. 

A  projection  on  the  moving  head  tripped  a  microswitch  (mounted 
on  one  leg  of  the  tripod  support  as  shown  in  Fig.  5.3) 
connected  to  an  alternate  on-off  type  relay,  which  starts  and 
stops  an  electric  interval  timer  (clock) .  The  relay  and  timer 
are  shown  on  top  of  the  machine  in  Fig.  5.4.  On  the  first  time 
around,  the  clock  was  started;  after  one  revolution  the  projec¬ 
tion  tripped  the  microswitch  again  and  the  clock  was  stopped. 

The  speed  measured  by  this  method  was  within  -  1.5%. 

(4)  Torque  Measurement  -  The  torque  was  measured 
with  a  "torque  cell"  connected  in  series  with  the  stationary 
head.  The  torque  cell  was  a  steel  spindle  with  two  sets  of 
strain  gages  on  it.  One  set  of  four  on  the  spindle  at 

45  degrees  to  the  axis  measured  the  shear  strain  during  torsion. 
The  other  set  of  four  on  the  spindle  parallel  to  the  axis  measured 
the  possible  strain  due  to  bending.  The  latter  set  was  installed 
for  preliminary  work  only;  it  was  used  to  check  for  bending 
during  the  test  and  to  check  the  effect  of  the  bending  on  the 
torque  readings.  It  was  found  that  the  bending  was  negligible 
and  did  not  influence  the  torque  measurements.  The  four  strain 
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Fig.  5.4  Another  View  of  the  Testing  Machine 


Fig.  5.5  Arrangement  for  Calibration  of  the  Torque  Cell 
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Fig,  5.6  Circuit  Diagram  for  Motor  Speed  Adjustments 


Fig.  5.7  Circuit  Diagram  for  Torque  Measurement 
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gages  on  the  spindle  at  45  degrees  to  the  axis  were  arranged 
so  that  two  were  in  compression  and  two  in  tension.  The  changes 
in  resistance  of  these  gages  were  measured  by  the  bridge-circuit 
shown  in  Fig.  5.7  and  the  bridge  was  balanced  by  adjusting  the 
decade  resistance  box  until  the  microvolt  meter  read  zero.  An 
amplification  of  the  input  voltage  was  taken  from  a  terminal  at 
the  back  of  the  microvolt  meter  and  fed  into  the  y-channel  of  an 
x-y  recorder  with  the  x-channel  in  time  sweep;  thus,  a  con¬ 
tinuous  record  of  the  torque  was  obtained  during  torsion  of 
each  specimen.  In  Fig.  5.2,  the  microvolt  meter,  the  6- volt 
power  supply,  and  the  decade  resistance  box  are  shown  on  the 
foreground  table,  and  the  x-y  recorder  is  on  the  background 
table.  The  torque  was  calibrated  by  actually  applying  increments 
of  pure  torque  on  the  stationary  head  in  the  manner  indicated 
in  Fig.  5.  5  and  recording  the  signals  on  the  x-y  recorder.  The 
length  of  the  arm  on  each  side  was  10.0  inches,  and  the  weights 
applied  on  each  side  ranged  from  2  lbs.  to  7  lbs.  in  increments 
of  one  or  lb. 

(5)  Temperature  Measurement  -  The  temperature  of  the 
specimen  during  test  was  measured  by  inserting  an  iron-constantan 
thermocouple  into  the  stationary  shoulder  of  the  specimen  as 
shown  in  Fig.  5.3,  and  the  temperature  was  measured  by  the  Leeds 
and  Northrup  potentiometer  shown  in  the  foreground  of  Fig.  5.2 
(with  cover  opened) .  The  temperature  was  not  monitored  continu¬ 
ously,  but  was  recorded  manually  at  fixed  intervals  only. 
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5.3  The  Test  Specimens 

The  test  specimens  were  one  inch  diameter  rods  with 
a  narrower  test  section  in  the  middle,  as  shown  in  Fig.  5.8. 

The  fillet  radii  at  both  ends  of  the  test  section  were 
approximately  3/32  inch.  The  "gage  length"  of  the  test  section 
was  then  defined  as  the  distance  between  the  two  points  of 
tangency  of  the  fillets.  This  definition  of  the  gage  length  is 
not  exactly  correct  because  the  twisting  of  the  material  does 
not  start  at  cross  sections  through  the  points  of  tangency  but 
at  cross  sections  through  the  fillets.  The  cross  sections  at 
which  twisting  starts  are  difficult  to  determine,  so  the  above 
definition  of  gage  length  was  taken  as  an  approximation.  The 
approximation  would  not  affect  the  results  of  the  experiment 
provided  the  definition  is  taken  consistently  throughout  the 
experiment;  this  was  done. 

The  test  specimens  were  made  from  one-inch  diameter 
rods  formed  from  commercially  pure  lead  billets  by  inverted 
extrusion.  It  is  usually  assumed  that  pure  lead  anneals  at  room 
temperature  as  quickly  as  it  workhardens.  However,  it  was  dis¬ 
covered  during  preliminary  tests  that  this  assumption  does  not 
appear  to  be  true  for  commercially  pure  lead  which  is  99.998% 
pure.  It  was  found  that  a  period  of  about  five  days  was  required 
for  the  extruded  rods  to  anneal,  and  that  higher  annealling 
temperatures  did  not  shorten  the  annealling  time.  The  diameter 
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Fig.  5.8  Drawing  of  a  Test  Specimen 
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Fig.  5.9  Typical  Test  Specimens 
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of  the  test  section  was  machined  to  within  -  0.001"  and  the 
taper  was  kept  to  within  0.0007"  per  inch. 

All  the  specimens  were  measured  with  a  Leitz  measur¬ 
ing  microscope  before  testing.  For  each  specimen,  the  diameters 
in  four  directions  and  at  three  points  along  the  test  section 
were  measured;  the  gage  length  was  measured  by  judging  the 
points  of  tangency  of  the  fillets  and  four  measurements  were 
made  around  the  circumference;  two  shoulder-to-shoulder  length 
measurements  were  made  around  the  circumference.  In  this  way, 
the  taper  of  the  test  section  was  checked,  the  average  diameters 
and  the  average  shoulder-to-shoulder  lengths  were  obtained  to 
within  ^  0.001",  and  the  average  gage  lengths  were  obtained  to 
within  -  0.002".  Two  typical  test  specimens  are  shown  in 
Fig.  5.9. 

About  130  5/8-inch  diameter  specimens  were  tested. 
About  70  of  these  specimens  were  used  for  preliminary  tests. 
Preliminary  tests  were  performed  to  determine  the  effect  of  gage 
length  on  the  torque  for  a  given  rate  of  twist  per  unit  length. 
It  was  found  that  change  of  gage  length  had  an  effect  on  the 
torque.  At  a  given  rate  of  twist  per  unit  length,  a  longer  gage 
length  gave  a  lower  torque.  This  is  probably  due  to  end  effects 
of  the  test  section.  The  material  near  the  non-deforming  shoul¬ 
ders  does  not  twist  as  much  as  the  material  away  from  the 
shoulders;  thus,  the  gage  length  is  affected.  Since  it  is 
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difficult  to  determine  quantitatively  this  end  effect,  to 
avoid  this  complexity,  a  constant  gage  length  of  approximately 
1.33"  (or  approximately  1.5"  shoulder-to-shoulder  length)  was 
used  only. 

Three  batches  of  commercially  pure  lead  were 
tested,  and  were  referred  to  as  Lead  A,  B,  and  C.  The 
suppliers  and  assays  are  given  below: 

(1)  Lead  A  -  Supplied  by  Canada  Metals  Ltd.  Its 
assay  obtained  by  approximate  analysis  is  0.02%  antimony, 

0.015%  arsenic,  and  by  difference  99.965%  lead. 

(2)  Lead  B  -  Supplied  by  Consolidated  Mining  and 
Smelting  Company  of  Canada,  Ltd.  The  supplier's  assay  is 
0.0002%  silver,  0.0002%  copper,  0.0003%  bismuth,  0.0002%  thallium, 
0.0001%  arsenic,  0.00002%  antimony  and  by  difference  99.998%  lead. 

(3)  Lead  C  -  Supplied  by  Consolidated  Mining  and 
Smelting  Company  of  Canada  Ltd.  The  supplier  claimed  that  this 
batch  is  the  same  as  Lead  B;  however,  slightly  different  results 
were  obtained  from  tests. 

One  hollow  specimen  with  3/4-inch  outer  test  diameter 
and  1/2-inch  inner  diameter  was  tested.  It  was  machined  from  a 
1-inch  diameter  lead  tube  which  was  formed  by  tube  extrusion. 

When  twisted,  it  started  to  buckle  and  a  "wringing"  phenomenon 
occurred.  The  same  phenomenon  was  observed  by  Swift  (24)  for 
hollow  specimens.  The  phenomenon  may  be  caused  by  instability  due 
to  slight  eccentricity  of  the  hole  or  slight  inhomogeneity  of 


the  material. 


,  Jy  i  ■  e  a.:  '  \i>  ...  n,;*  b  \el  ol  .  Luo  »:  >  . 

•  \  • 

4.  /  '  1  1  •;  -*.r;  io  ;  ■’  ; 

«  ;  .  f-  .  -  •  • 

WQlg  ;  tfe  ip  >-  »  i  XO.  •  *  i£ 

1.  -■  '.  if»-.  -  •  t 

4 .»  .:.4  J.6r -4.:-0' ,  ■  _  .  [C  ..-9  st  do  2J.S 

b  '  .  ;■  •  r;"^ v  •  .n:  ■*.’  -1 

!  ;  .  •  1  \  1  a  -  t  or.  . 

>  0,  '  ,  t  ) 00.  >qq:  ►  *?SOOO.O  ,  .'SCOC  (• 

0  is'* ■  •  •  '  1  .ir  cr;  ■  .000  0  ,  .n  -i?xb  ^XOO.). 

;  r\.  .  fc0V*->&-  -»  u:W  i  .  •  JOflOf .  ..  '■/’  “  ..  ■  '  U  DDO 

f€  ::  \f  \  •  if'  ...  ■  ..  ..  :•  -  X  i 


76 


5.4  Experimental  Procedure 

The  experiments  consisted  of  twisting  the  speci¬ 
mens  at  various  rates  of  twist  per  unit  gage  length  over  a 
range  of  0.75  rpm/inch  to  11  rpm/inch. 

The  experimental  procedure  for  testing  each  specimen 
started  with  the  calibration  of  the  torque  cell  by  first  clamping 
the  calibration  arm  on  the  stationary  head  as  shown  in  Fig.  5.  5.. 
Before  weights  were  applied  on  the  arm,  the  bridge  circuit  was 
first  balanced,  the  sweep  of  the  torque  measuring  x-y  recorder 
was  started,  and  the  pen  was  zeroed.  Then  increments  of  weight 
were  applied  to  give  torques  of  70,  80,  and  90  inch-pounds. 

After  each  increment  of  weight  was  applied,  a  slight  pressure 
was  applied  and  released  manually  on  the  arm  in  the  direction 
that  increases  the  torque  on  the  head.  This  reversed  the 
direction  of  the  impending  motion  of  the  pulleys  and  consequently 
the  direction  of  the  frictional  forces.  A  slight  pressure  was 
applied  and  released  again  on  the  arm  in  the  opposite  direction 
returning  the  pulleys  to  the  original  impending  motion.  The 
average  of  the  traces  obtained  from  the  two  impending  motions 
was  taken  as  the  true  torque  calibration.  However,  friction  in 
the  pulleys  was  not  significant  and  negligible  difference  between 
the  two  traces  was  noticed.  A  typical  calibration  trace  is  shown 
in  Fig.  5.10.  Calibration  of  the  torque  cell  was  done  for  each 
specimen  tested  in  order  to  eliminate  the  error  of  possible 
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Ficr.  5.10  Typical  Toraue  and  Calibration  Curves 
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change  of  calibration  between  tests;  also  it  was  much  more  con¬ 
venient  to  have  the  calibration  on  the  same  chart  sheet  as  the 
torque  curve.  After  the  calibration,  the  machine  was  started 
and  the  speed  was  adjusted  to  the  required  deflection  on  the 
speed  monitoring  x-y  recorder;  the  required  deflection  to  obtain 
a  given  rate  of  twist  during  the  test  was  determined  by  prelimi¬ 
nary  tests.  After  the  speed  adjustment  was  made,  a  "no  load" 
trace  was  made  on  the  x-y  recorder  to  check  the  speed  adjustment. 
Then  the  machine  was  stopped,  the  calibration  arm  was  removed 
and  replaced  by  the  specimen,  and  the  thermocouple  was  inserted 
in  the  stationary  shoulder  of  the  specimen  as  near  as  possible 
to  the  test  section.  Before  the  lower  shoulder  of  the  test 
specimen  was  clamped  in  the  moving  head,  the  balance  of  the 
bridge  circuit  was  checked,  the  torque  measuring  x-y  recorder 
was  zeroed  and  started  to  sweep,  and  the  initial  temperature 
of  the  specimen  was  recorded  from  the  potentiometer.  Then  the 
specimen  was  clamped  in  the  moving  head  making  certain  that 
movement  clearance  was  provided  in  the  slider  at  the  coupling  of 
the  moving  head  shaft  and  that  the  specimen  is  not  pre-loaded, 
that  is,  the  torque  measuring  x-y  recorder  still  reads  zero. 

With  both  x-y  recorders  reading  zero  and  sweeping,  the  machine 
was  started.  Then  the  rate  of  twist  curve  and  torque  curves 
were  simultaneously  recorded  on  the  x-y  recorders.  Also  simul¬ 
taneously,  the  time  for  every  second  revolution  of  the  moving 
head  was  obtained  from  the  clock.  The  time  intervals  were 
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recorded  manually.  As  each  time  the  relay  coil  was  energized 
to  start  and  stop  the  clock,  the  magnetic  induction  from  the 
coil  caused  a  blip  on  the  torque  curve;  these  blips  conveniently 
indicated  the  number  of  revolutions  but  did  not  affect  the 
accuracy  of  the  torque  measurement.  At  every  click  of  the 
relay,  the  temperature  was  recorded  manually,  giving  the  tem¬ 
perature  of  the  specimen  at  every  revolution.  (A  typical 
torque  trace  is  shown  in  Fig.  5.10) .  While  the  specimen  was 
being  twisted,  it  was  checked  constantly  for  transverse  shear 
failure.  If  this  occurred,  the  machine  was  stopped.  If  no 
shear  failure  occurred,  the  machine  was  stopped  after  1~  to  2 
minutes  (or  4  to  20  revolutions)  of  twisting,  depending  on  the 
speed  of  twist.  After  the  machine  was  stopped,  the  moving  head 
jaws  were  loosened  freeing  the  torque  cell  from  torque;  then  the 
torque  trace  was  checked  to  see  if  the  pen  returned  to  zero. 
After  the  specimens  were  tested,  the  shoulder-to-shoulder 
lengths  and  the  diameters  were  measured  again  by  the  measuring 
microscope . 

5.5  Experimental  Results  and  Discussion 

All  the  torque  curves  show  a  bump  at  the  beginning  of 
the  test,  then  the  torque  drops  to  a  minimum  and  increases 
slightly  during  the  rest  of  the  test.  No  explanation  is 
offered  for  this  initial  bump.  Such  bumps  were  obtained  by 
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Loizou  and  Sims  (20)  and  by  Kobayashi,  et  al  (21)  in  compression 
tests.  The  slight  subsequent  increase  of  torque  during  the 
test  is  due  to  the  increase  in  diameter  of  the  specimen,  which 
results  from  the  shortening  of  the  specimen.  The  torque  used 
for  correlation  with  the  theoretical  result  (5.6)  was  taken  as 
the  minimum  torque  after  the  bump,  that  is  before  any  signifi¬ 
cant  shortening  and  increase  in  diameter  occurred. 

The  rate  of  twist  was  calculated  from  the  clock 
reading  for  the  revolution  in  which  the  minimum  torque  occurred. 
The  speed  monitoring  curve  was  checked  to  ensure  that  the  speed 
during  that  revolution  was  uniform.  The  number  of  revolutions 
for  each  test  was  obtained  by  counting  the  blips  on  the  torque 
curve,  and  the  change  of  length  of  each  specimen  was  obtained 
from  the  difference  between  the  shoulder-to-shoulder  lengths 
before  test  and  after  test.  (Continue  on  page  87). 
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Table  5.1 

Summary  of  Data  for  Lead  A 


Gage 

Length 

Rate  of 
Twist 

CD 

L 

Torque 

Change  in  Temp. 

(°F) 

L 

(in. ) 

CD 

(rpm) 

(rpm/in. ) 

(in. /lbs. ) 

At  Min . 

Torque 

At  End 
of  test 

1.353 

2.02 

1.49 

77.8 

1.9 

2.5 

1.362 

2.06 

1.51 

78.3 

2.7 

2.8 

1.340 

2.67 

1.99 

78.8 

3.5 

4.4 

1.309 

3.57 

2.73 

81.5 

1.9 

5.2 

1.362 

4.55 

3.34 

82.2 

5.4 

6.9 

1.332 

5.17 

3.88 

82.5 

5.6 

8.0 

1.378 

5.35 

3.88 

83.1 

6.3 

8.5 

1.309 

6.00 

4.59 

85.0 

— 

— 

1.338 

7.23 

5.40 

87.0 

9.7 

12.3 

1.345 

7.31 

5.44 

89.5 

11.8 

14.3 

1.327 

7.70 

5.80 

86.3 

9.6 

11  o  8 

1.333 

8.11 

6.09 

86.9 

10.5 

13.7 

1.336 

8.82 

6.60 

88.6 

12.1 

15.8 

1.336 

8.95 

6.70 

87.0 

10.7 

15.5 

1.331 

9.23 

6.94 

88.3 

15.3 

18.0 

1.316 

9.53 

7.25 

86.5 

14.0 

17.6 
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Table  5.2 

Continuation  of  Summary  of  Data  for  Lead  A 


Speci¬ 
men  No. 

Gage 
Length 
Before 
Test 
(in. ) 

Shoulder-to-Shoulder 
Length  (in.) 

AL 

L 

(%) 

No .  of 
Revs . 
Twisted 

ND 

L 

Before 

Test 

After 

Test 

54 

1.353 

1.5081 

1.4943 

1.02 

4.60 

2.12 

52 

1.362 

1.5083 

1.4922 

1.18 

4.70 

2.16 

53 

1.340 

1.5020 

1.4811 

1.56 

8.40 

3.92 

56 

1.378 

1.5200 

1.4977 

1.62 

8.70 

3.95 

55 

1.362 

1.5108 

1.4883 

1.65 

9.60 

4.41 

4 

1.309 

1.5037 

1.4832 

1.57 

12.35 

5.90 

34 

1.338 

1.5093 

1.4817 

2.06 

16.50 

7.71 

5 

1.309 

1.5003 

1.4672 

2.53 

16.95 

8.09 

29 

1.345 

1.5122 

1.4752 

2.75 

17.20 

7.99 

33 

1.327 

1.5035 

1.4642 

2.96 

20.85 

9.82 

31 

1.336 

1.5088 

1.4713 

2.81 

21.60 

10.10 

30 

1.331 

1.5059 

1.4626 

3.25 

23.10 

10.85 

23 

1.336 

1.5103 

1.4656 

3.35 

24.05 

11.25 

22 

1.316 

1.4999 

1.4534 

3.53 

26.70 

12.68 

*• 

- toeqa 

J-  •- 

. 

■ 

T 

fT  £ 

i>2.1 

:  {  ' 

:•.  J 

'  • 

Sd.I 

t38*. C 

•Oi-.  • 

£>  ■ 

t?oa.i 

••  . . 

r Z.XL 

• 

is. s: 

£IV^.X 

' 

0 

cao  a.x 

• 

eo.K 

dfcdKJ 

eoxa.i 

£S 

72 

71 

75 
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76 
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88 
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96 
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Table  5.3 

Summary  of  Data  for  Lead  B 


Gage 

Length 

L 

(in.) 

Rate  of 
Twist 

0) 

(rpm) 

0) 

L 

(rpm/in. ) 

Torque 

(in. /lbs . ) 

Change 

( 

At  Min. 
Torque 

in  Temp. 

°F) 

At  End 
of  test 

1.331 

2.00 

1.50 

79.0 

2.0 

2.9 

1.345 

2.03 

1.51 

78.0 

2.5 

2.9 

1.341 

2.64 

1.97 

78.4 

2.8 

3.7 

1.331 

3.61 

2.71 

81.7 

3.4 

5.1 

1.352 

4.38 

3.24 

81.2 

3.9 

6.8 

1.333 

5.05 

3.80 

82.6 

5.6 

8.3 

1.362 

6.00 

4.40 

84.6 

6.7 

10.7 

1.381 

6.19 

4.48 

85.0 

6.5 

9.7 

1.337 

6.67 

4.99 

84.5 

8.0 

11.1 

1.335 

6.67 

5.00 

83.9 

6.2 

10.8 

1.335 

6.82 

5.11 

85.5 

7.2 

10.8 

1.324 

7.06 

5.33 

86.0 

8.1 

12.3 

1.359 

8.00 

5.89 

85.8 

9.6 

11.8 

1.347 

8.82 

6.55 

83.9 

8.2 

12.5 

1.343 

10.00 

7.44 

84.7 

9.6 

17.7 

1.346 

11.32 

8.42 

86.0 

8.8 

18.0 

1.362 

12.50 

9.19 

86.0 

12.3 

22.5 

1.331 

12.24 

9.20 

86.3 

9.3 

15.6 

1.365 

13.64 

10.00 

85.8 

11.3 

18.7 

1.326 

15.00 

10.92 

86.0 

12.3 

21.0 

1.373 

15.40 

11.20 

87.3 

12.7 

25.7 

-  .. 
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Table  5.4 


Continuation  of  Summary  of  Data  for  Lead  B 


Speci¬ 
men  No. 

Gage 

Length 

Before 

Test 

(in.) 

■  . . L 

Shoulder-to-Shoulder 
Length  ( in . ) 

AL 

L 

(%) 

No .  of 
Revs . 
Twisted 

ND 

Before 

Test 

After 

Test 

L 

72 

1.331 

1.4966 

1.4820 

1.10 

4.45 

2.09 

71 

1.345 

1.5009 

1.4884 

0.93 

4.65 

2.16 

75 

1.341 

1.5041 

1.4851 

1.42 

6.15 

2.87 

85 

1.331 

1.5016 

1.4834 

1.37 

6.25 

2.93 

89 

1.358 

1.5203 

1.5037 

1.22 

6.40 

2.95 

76 

1.352 

1.5132 

1.4915 

1.60 

8.85 

4.09 

88 

1.381 

1.5476 

1.5232 

1.74 

10.00 

4.53 

77 

1.333 

1.4962 

1.4650 

2.34 

12.35 

5.79 

86 

1.337 

1.5020 

1.4700 

2.39 

12.80 

5.98 

78 

1.362 

1.5237 

1.4871 

2.69 

13.25 

6.08 

90 

1.324 

1.4882 

1.4479 

3.04 

13.60 

6.42 

92 

1.359 

1.5146 

1.4746 

2.94 

14.45 

6.65 

79 

1.335 

1.5060 

1.4652 

3.06 

14.70 

6.88 
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Table  5 . 5 


Summary  of  Data  for  Lead  C 


Speci¬ 
men  No. 

Gage 

Length 

Rate  of 
Twist 

CD 

L 

Torque 

Change  in  Temp. 
(°F) 

L 

(in. ) 

CD 

(rpm) 

(rpm/in.) 

( in ./lbs . ) 

At  Min. 
Torque 

At  End 
of  test 

102 

1.346 

0.985 

0.73 

74.0 

1.6 

1.6 

106 

1.352 

1.004 

0.74 

75.3 

1.5 

1.5 

103 

1.329 

1.95 

1.47 

76.8 

2.7 

3.0 

104 

1.321 

1.97 

1.49 

78.0 

2.4 

3.0 

111 

1.353 

3.06 

2.26 

82.5 

4.0 

4.7 

110 

1.372 

3.11 

2.27 

80.5 

4.2 

5.5 

113 

1. 313 

4.14 

3.15 

83.0 

4.8 

6.5 

112 

1.337 

4.29 

3.20 

83.2 

5.0 

6.8 

119 

1.355 

5.09 

3.75 

85.7 

5.5 

7.5 

114 

1.348 

5.13 

3.80 

83.8 

5.5 

8.8 

120 

1.346 

5.94 

4.41 

87.2 

7.1 

9.3 

124 

1.332 

6.12 

4.60 

86.1 

6.9 

9.3 

125 

1.367 

6.98 

5.10 

86.4 

5.7 

11.7 

118 

1.345 

7.32 

5.45 

85.8 

— 

— 

126 

1.385 

7.90 

5.70 

87.0 

7.3 

11.0 

100 

1.340 

8.95 

6.69 

88.2 

9.8 

14.3 

117 

1.353 

9.38 

6.93 

87.6 

8.7 

11.4 

107 

1.335 

10.34 

7.75 

87.2 

8.7 

11.6 

108 

1.329 

10.53 

7.93 

88.5 

9.1 

14.5 

101 

1.337 

11.54 

8.64 

89.2 

11.7 

17.5 

121 

1.356 

12.00 

8.85 

88.9 

12.3 

17.3 

115 

1.363 

13.33 

9.78 

89.4 

9.9 

15.1 

116 

1.332 

13.05 

9.80 

89.1 

9.6 

16.0 

123 

1.377 

15.00 

10.90 

89.4 

— 

— 

122 

1.358 

15.00 

11.03 

87.6 

12.0 

19.  .6 

. 
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Table  5.6 

Continuation  of  Summary  of  Data  for  Lead  C 


Speci¬ 
men  No. 

Gage 
Length 
Before 
Test 
(in. ) 

Shoulder-to 

Length 

-Shoulder 

(in.) 

AL 

L 

(%) 

No .  of 
Revs . 
Twisted 

ND 

Before 

Test 

After 

Test 

L 

102 

1.346 

1.5108 

1.4980 

0.95 

3.5 

1.63 

106 

1.352 

1.5076 

1.4934 

1.05 

4.0 

1.85 

104 

1.321 

1.4980 

1.4847 

1.01 

4.4 

2.08 

103 

1.329 

1.5085 

1.4925 

1.20 

4.5 

2.12 

110 

1.372 

1.5145 

1.4942 

1.48 

6.2 

2.82 

113 

1.313 

1.4826 

1.4660 

1.26 

6.3 

3.00 

111 

1.353 

1.5018 

1.4863 

1.15 

6.7 

3.10 

119 

1.355 

1.5094 

1.4824 

1.99 

6.9 

3.18 

117 

1.353 

1.5084 

1.4740 

2.54 

7.1 

3.28 

107 

1.335 

1.4951 

1.4673 

2.08 

7.2 

3.37 

120 

1.346 

1.4998 

1.4615 

2.85 

8.1 

3.76 

124 

1.332 

1.5044 

1.4516 

3.96 

8.6 

4.04 

109 

1.335 

1.5155 

1.4742 

3.09 

8.7 

4.07 

118 

1.345 

1.5027 

1.4627 

2.97 

8.9 

4.14 

108 

1.329 

1.4959 

1.4579 

2.86 

8.9 

4.19 

115 

1.363 

1.5262 

1.4651 

4.48 

9.6 

4.40 

121 

1.356 

1.5141 

1.4468 

4.96 

10.7 

4.93 

125 

1.367 

1.5257 

1.4600 

4.80 

10.9 

4.98 

116 

1.332 

1.4949 

1.4260 

5.17 

10.8 

5.07 

100 

1.340 

1.5082 

1.4525 

4.16 

11.1 

5.18 

122 

1.358 

1.5136 

1.4478 

4.85 

11.5 

5.29 

101 

1.337 

1.5089 

1.4398 

5.17 

12.7 

5.94 

•  • 
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Table  5.7 


Estimated  Accuracy  of  Data 


Data 


Accuracy 


Gage  Length 
Rate  of  Twist 
Torque 
No.  of  Revs. 
Change  of  Length 


t  0 . 002 " 

*  1.596 

i  0.3  in- lbs 
±  0.5% 

-  0.003" 


For  specimens  that  were  twisted  at  rates  of  twist 
per  unit  gage  length  from  approximately  5.0  rpm/inch  to  5.5  rpm/inch, 
the  torque  curves  were  wavy,  as  shown  in  Fig. 5.11.  The  possibility 
that  this  was  due  to  a  machine  characteristic  was  investigated 
and  eliminated,  and  no  satisfactory  explanation  could  be  found 
for  this  irregularity.  Accurate  torques  could  not  be  obtained 
for  this  range  from  the  torque  curves  because  of  the  waviness,  so 
only  the  average  torques  were  taken  as  approximations  in  the  range 
5.0  rpm/inch  to  5.5  rpm/inch. 

All  the  data  obtained  for  each  specimen  are  pre¬ 
sented  in  Tables  5.1  to  5.6.  The  estimated  accuracy  of  the  data 
are  given  in  Table  5.7.  The  curves  plotting  torques  against  rates 
of  twist  per  unit  length  (-r)  and  plotting  percentage  shortening 
of  the  original  gage  length  against  the  dimensionless  parameter 
”,  where  D  is  the  diameter  of  the  specimen  and  L  is  the 
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Fig.  5.11  A  Typical  Wavy  Torque  Curve 
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original  gage  length,  are  shown  in  Figs.  5.12  to  5.17  for  the 
three  batches  of  pure  lead. 

The  torques  appear  to  have  a  linear  relation  with  rates 
of  twist  per  unit  length  for  ~<5.0  rpra/inch  ,  and  appear  to  be 
independent  of  ~  for  —>5.5  rpm/inch.  In  the  region  between 
5.0  rpm/inch  and  5.5  rpm/inch,  the  data  are  quite  scattered; 
consequently,  no  conclusion  can  be  made.  There  is  a  possibility 
that  a  critical  region  exists  near  the  point  where  the  rate 
dependent  characteristics  of  lead  change.  This  region  corresponds 
approximately  to  a  maximum  shear  strain  rate  (engineering  defini¬ 
tion)  in  the  specimen  of  about  0.16  sec  ^ .  In  the  investigation 
of  Loizou  and  Sims  (20)  using  the  compression  test,  no  such  criti¬ 
cal  region  existed.  This  critical  phenomenon  may  be  inherent  to  the 
torsion  test,  which  subjects  the  material  to  shear  only. 

Straight  lines  were  fitted  by  the  method  of  least 
squares  to  the  torque  versus  ^  curves  for  the  range  of 
~  <  5.0  rpm/inch.  For  each  lead,  equation  (5.16)  in  the  form 

Li 

T  =  A  +  An  (7^)  was  obtained.  The  constants  A  and  An  ,  are 

O  L  Lj  OX 

shown  in  Table  5.8.  From  these  constants,  the  coefficients  of 
viscosity  t)  and  yield  stresses  in  shear  k  were  calculated,  and 
are  also  shown  in  Table  5.8.  An  estimate  of  the  error  of  r\  may 
be  obtained  by  an  estimate  of  the  error  of  A^ ,  which  is  the  slope 
of  the  fitted  straight  line.  From  a  statistics  reference  (26) , 
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Fig.  5.12  Torque  vs.  Rate  of  twist  per  Unit  Length  for  Lead 
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Fig.  5  13  Torque  vs.  Rate  of  Twist  per  Unit  Length  for  Lead 
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.14  Torque  vs.  Rate  of  Twist  per  Unit  Length  for  Lead  C 
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Fig.  5.15  Percentage  Contraction  of  Gage  Length  vs.  Dimensionless 

Parameter  ND/L  for  Lead  A 
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Fig.  5.16  Percentage  Contraction  of  Gage  Length  vs.  Dimensionless 

Parameter  ND/L  for  Lead  3 
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Fig.  5.17  Percentage  Contraction  of  Gage  Length  vs.  Dimensionles 

Parameter  ND/L  for  Lead  C 
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the  standard  errors  of  the  slopes  were  calculated.  Since  the 
error  of  the  radius  of  the  specimen  is  negligible,  the  percentage 
error  of  slope  would  then  be  the  percentage  error  of  q.  (Refer 
to  equation  (5.6)).  The  estimated  percentage  errors  of  q  for 
the  three  batches  of  lead  are  shown  also  in  Table  5.8. 


Table  5.8 


Coefficients  of 

Viscosity 

and  Yield 

Stresses 

in  Shear 

A 

0 

Ai 

T) 

(psi-sec) 

k 

...  (psi) 

Estimated 
Error  of 

% 

T) 

Lead  A 

74.85 

2.152 

1372 

1170 

7% 

Lead  B 

75.19 

2.047 

1305 

1176 

13% 

Lead  C 

73.02 

3.135 

1998 

1142 

8% 

It  can  be  seen  from  Table  5.8  that  the  rate  sensitiv¬ 
ity  of  commercially  pure  lead  varies  greatly  with  small  changes  of 
composition.  This  contributed  to  the  difficulties  of  the  experi¬ 
ments.  The  variation  of  the  rate  sensitivity  may  be  due  to  the 
extent  of  annealling  of  the  lead,  the  casting  temperature,  or  the 
rate  of  cooling  when  the  lead  was  cast. 

The  data  given  in  Figs.  5.15  and  5.16  show  a  linear 
correlation  between  the  percentage  shortening  and  the  parameter 
~  9  so  straight  lines  were  fitted  to  the  data.  In  Fig.  5.17,  a 
best  fit  straight  line  was  drawn  for  the  data,  but  the  scattering  is 
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such  that  no  definite  conclusion  can  be  made.  Irregularity  due  to 
friction  at  the  slider  coupling  of  the  shaft  of  the  moving  head  is 
a  possible  reason  for  the  scattering  of  the  data.  Swift  (24) 
twisted  one  lead  specimen;  the  curve  plotting  percentage  shortening 
against  shear  strain  at  the  outer  diameter  during  torsion 
showed  also  a  linear  relationship. 

The  slopes  of  the  straight  lines  fitted  to  data  in 
Figs.  5.15  to  5.17  are  respectively  0.237,  0.424,  and  1.204. 

This  variation  of  the  slopes  shows  that  the  shortening  appears  to 
be  quite  sensitive  to  the  metallurgy  of  the  pure  leads. 

Swift  (24)  gave  an  explanation  for  the  change  of  length 
of  material  during  pure  torsion  by  qualitative  considerations  of 
the  differential  grain  orientations.  Hill  (25)  suggested  that 
the  change  of  length  is  due  to  both  anisotropy  and  differential 
orientations  of  the  grains,  and  derived  an  equation  for  the  change 
of  length  assuming  the  change  is  solely  due  to  anisotropy. 

5 . 6  Cone lusions  and  Recommendations 

The  results  of  this  experimental  investigation  show  that 
the  torsion  test  can  be  used  to  assess  strain  rate  sensitivity  of 
materials,  although  the  change  in  length  of  the  specimen  and  the 
bump  in  the  torque  curve  may  influence  the  assessment.  The 
linear  dependence  of  the  torque  on  the  rate  of  twist  per  unit  length 
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shows  that  commercially  pure  lead  behaves  according  to  the 
Bingham  visco-plastic  model;  however,  this  is  true  only  under  a 
certain  critical  region;  beyond  the  critical  region,  the  torque  is 
independent  of  rate  of  twist  per  unit  length.  The  results  show 
that  reasonable  coefficients  of  viscosity  can  be  obtained  by 
the  torsion  test. 

The  critical  region  is  of  interest;  however,  no  general 
conclusion  can  be  made  and  further  work  is  needed  to  investigate 
this  phenomenon.  Further  work  may  also  be  done  by  varying  the 
radius  of  the  specimen  to  verify  further  that  pure  lead  behaves 
according  to  the  Bingham  visco-plastic  model.  For  future  work, 
better  metallurgical  control  of  the  pure  lead  is  suggested. 

The  correlation  between  the  shortening  of  the 
specimens  and  the  parameter  proportional  to  the  number  of  revolutions 
twisted  shows  a  tendency  to  be  linear;  however,  because  of  the 
scattering  of  the  data,  no  definite  conclusion  can  be  made. 

Further  work  may  be  done  by  improving  the  method  of  allowing  the 
specimen  to  change  length  during  twist  and  the  method  of  measur¬ 
ing  the  change  of  length,  for  example,  through  continuous  measure¬ 
ment  by  electrical  displacement  transducers. 
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APPENDIX 

COMPUTER  PROGRAMS  AND  CALCULATED  RESULTS 


The  computer  programs  are  written  in  Fortran  IV  Language 
for  an  IBM  7040  Digital  Computer „  The  following  nomenclatures 
are  used  throughout  all  the  computer  programs  and  symbols 
not  listed  are  particular  to  each  programs 


E 

H 

N 

P(I,J) 

SI 

sx 

SY 

ST 

T(I,J) 

TQQ 


grid  width 

one  plus  the  number  of  grid  divisions 

I 

s 


(i,  j) 


xz 


=  s 


yz 

T  (for  Program  B) 
t  (except  Program  B) 
T 


:■  .f  •  ■  V  ■ 

’ 

'o/i.  <■'>  '  lz  '<  (?  c  “  ■  -17 


I  “1  r  ' 


Fic 
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COMPUTER  PROGRAM  a 

#♦ 

DIMENSION  P ( 51  *  51  )  *  T ( 51 >  51  ) *TDX( 51*51) >TDY(51*51) 

READ  (5*100)  E  *  H  >  N 
NN=N-1 
HH  =  H*H 
H2  =  HH/ 2 • 

H6=H/6. 

DO  5  I  =  1  *  N 
5  P(  I  *  N ) =0*  0 
Y=  1  •  0 
J=  N+l 
10  J= J-l 
K  =  J-l 

T  (  J  *  J  )  =  1  •  0 
TDX( J, J)=-2.*E 
TD Y ( J ♦ J ) =  2  •  *E 
X=-H 

DO  20  1=1 *K 
X=  X  +  H 

T (  I  >J)=1.+2.*E*(Y-X) 

TD X ( I  * J)  =  -2.*E 
TD Y ( I  * J  )=  2  •  *E 
TA  =  SIN(P(I,J)) /COS(P (  I  ,  J)  ) 

20  P(  I  * J-l )  =  P( I  *  J ) +TA* ( P ( 1+1 , J)-P( I , J)  ) -H* ( TDX ( I » J )-TA*TDY( I , J) ) /T(  I  * 
1  J) 

Y=  Y-H 

IF  (J.GT.2)  GO  TO  10 
WRITE  (6*101) 

TQ= ( l.+E*H/2. )*HH*H/32. 

Y=  0  • 

DO  45  J  =  2  * NN 
Y=  Y+H 
X=  0  , 

1=1 

SX  =  T ( I  * J)*COS(P( I  * J ) ) 

SY  =  T( I , J ) *5 1 N ( P ( I ,J)  ) 

TQ=TQ+(SX*Y+SY*H/4. )*H2  ^ 

WRITE  (6*102)  X,Y,P( I ,J) »SX*SY*T( I*J) 

25  IF  (J.EQ.2)  GO  TO  42 
K=  J-l 

DO  40  I =2 , K 
X=X+H 

SX  =  T ( I , J)*COS(P( I  *  J )  ) 

5  Y  =  T ( I » J ) *5  I N { P ( I , J) ) 

TQ=TQ+ ( SX*Y+SY*X ) *HH 
40  WRITE  (6*102)  X  * Y , P ( I , J ) * SX » SY , T ( I  * J ) 

42  I  =  J 
X=  X+H 

SX  =  T ( I  ,  J)*COS(P(  I , J)  ) 

S Y=  T ( I  *  J ) *S I N ( P ( I ,J)  ) 

SXU=SX*5. /6.+T( I*J+1 )*COS(P( I*J+1) )/6. 

TQ=TQ+(SXU*(Y+H6)+SY*(X-H6) )*H2 
45  WRITE  (6*102)  X » Y * P ( I , J ) * SX , SY , T ( I , J ) 

J  =  N 


H  I 


<X-Y>  *3**$+.r*(  tl  1 


Ml  .  5  >*  -  "0  \(  . «•  I  »-•  M< 


0  0T  CD  <S.T0* U  ! 


.SF\H*HM*  ?  ,  S  VH v  M  .1  1  ■-’  r 
■  ’  .  MM 


•  ( L ,  I  «  *  {  ,  T 

SH*  {  .>\H»Y<  «-Y *X5>  u?T «0T 

(  <U  :  ><UeO}*Ct«  )T.X? 


un  «<<*U  M^YtX  tSOX,^)  -ir!*V  0 


(  L«  )<M  2CO*<L. I ; T^X? 

(  tl,  i  )c  H-  !«  Ml  .  I  )  -  " 

,  i  1  )  -v  >*  UI>T  +  ^  .  '  --  • 

(  ~X)*Y  Md*  ’  Y  -  ))  ‘  )  - 

• :  >t< v’,  .a.n^«YtX  t$o.r«d)  tiji  * 
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Y=  1  •  0 

YY  =  1 • -H/4 • 

X=-H 

00  50  1=1, NN 
X=  X+H 
SX=T ( I ,  J) 

SY  =  0  • 

A  A  =  H2 

IF  (  I.EQ. 1 )  AA  =  H2 / 2 • 

TQ=TQ+SX*AA*YY 

50  WRITE  (6*102)  X  * Y , P (  I  ,  J ) , SX * SY , T (  I , J ) 
TQ=TQ+HH*YY/8  • 

TQQ=TQ*8. 

WRITE  (6*103)  TQ^ 

WRITE  (6*104)  TO 


80  STOP 

TOO  FORMAT  ( F 1 0 , 3 , F9 . 3  *  I  6  ) 
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FORMAT 

(75H 

I  J  X 

Y  P(I*J) 

1  SX 

SY/  ) 

102 

FORMAT 

( 1 1 X  *  2  F8 

•3*4F13.8) 

103 

FORMAT 

(31HK 

DIMENSIONLESS 

TORQUE  = , F 1 2 • 8 ) 

104 

FORMAT 

END 

(  3 1 H 

ONE-EIGHTH  OF 

IT  =  *  F12 • 8  ) 

DATA  CARD  - 


» J  ) 


0.100 


0.020 


51 


.f  V  It 


I-./ 

• !  \  h  -  ■  {  .  )■  .i )  -r 

-*■  n)  r, 

<i  ■  w  -  t.  -  - 

■-  ..  -1  M  I  i 
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RESULTS  FROM  PROGRAM  A 


X 

y 

4) 

-s 

xz 

Syz 

T 

0.000 

0.040 

0.19189749 

0.98949722 

0.19224766 

1.00799997 

0.040 

0.040 

0.19346200 

0.98134452 

0.19225745 

1 . 000000^0 

0.000 

0.080 

0.18225625 

0.99917221 

0.18414889 

1.01599997 

0.040 

0.080 

0.18372708 

0.99103492 

0.18415675 

1.00799997 

0.080 

0.080 

0.18522210 

0.98289536 

0.18416484 

1.00000000 

0.000 

0.120 

0.17279582 

1.00875045 

0.17606368 

1.02399995 

0.040 

0.120 

0.17417647 

1.00062745 

0.17606987 

1.01599997 

0.080 

0.120 

0.17557956 

0.99250243 

0.17607624 

1 .00799997 

0.120 

0.120 

0.17700567 

0.98437535 

0.17608282 

1 • nooo^o^^ 

O.'OOO 

0.160 

0.16350959 

1.01823520 

0.16799100 

1.03199996 

0.040 

0.160 

0.16480335 

1.01012541 

0.16799575 

1.02399997 

0.080 

0.160 

0.16611791 

1.00201383 

0.16800063 

1.01599997 

0.120 

0.160 

0.16745380 

0.99390040 

0.16800568 

1.00799997 

0.160 

0.160 

0.16881154 

0.98578513 

0.16801090 

1 . oooooooo 

0.000 

0.200 

0.15439133 

1.02762947 

0.15992984 

1.03999996 

0.040 

0.200 

0.15560132 

1.01953186 

0.15993335 

1.03199996 

0.080 

0.200 

0.15683055 

1.01143269 

0.15993696 

1.02399997 

0.120 

0.200 

0.15807952 

1.00333188 

0.15994071 

1.01599997 

0.160 

0.200 

0.15934869 

0.99522945 

0.15994458 

1.00799997 

0.200 

0.200 

0.16063857 

0.98712534 

0.15994858 

1 • oooonono 

0.000 

0.240 

0.14543512 

1.03693614 

0.15187926 

1.04799996 

0.040 

0.240 

0.14656428 

1.02884976 

0.15188172 

1.03999996 

0.080 

0.240 

0.14771123 

1.02076201 

0.15188426 

1.03199996 

0.120 

0.240 

0.14887639 

1.01267285 

0.15188689 

1.02399997 

0.160 

0.240 

0.15006021 

1.00458222 

0.15188963 

1.01599997 

0.200 

0.240 

0.15126314 

0.99649014 

0.15189246 

1.00799997 

0.240 

0.240 

0.15248564 

0.98839657 

0.15189539 

1 • oooooooo 

0.000 

0.280 

0.13663535 

1.04615792 

0.14383839 

1.05599996 

0.040 

0.280 

0.13768650 

1.03808186 

0.14383996 

1 . 04799996 

0.080 

0.280 

0.13875403 

1.03000460 

0.14384159 

1.03999996 

0.120 

0.280 

0.13983835 

1.02192613 

0.14384330 

1.03199997 

0.160 

0.280 

0.14093985 

1.01384639 

0.14384507 

1.02399997 

0.200 

0.280 

0.14205895 

1.00576537 

0.14384692 

1.01599997 

0.240 

0.280 

0.14319607 

0.99768303 

0.14384885 

1.00799997 

0.280 

0.280 

0.14435165 

0.98959937 

0.14385086 

1.00000000 

0.000 

0.320 

0.12798673 

1.05529735 

0.13580640 

1.06399996 

0.040 

0.320 

0.12896251 

1.04723078 

0.13580723 

1.05599996 

0.080 

0.320 

0.12995336 

1.03916314 

0.13580811 

1.04799996 

0.120 

0.32C 

0.13095964 

1.03109447 

0.13580904 

1.03999997 

0.160 

0.320 

0.13198170 

1.02302470 

0.13581002 

1.03199997 

0.200 

0.320 

0.13301992 

1.01495383 

0.13581105 

1.02399997 

0.240 

0.320 

0.13407470 

1.00688181 

0.13581214 

1.01599997 

0.280 

0.320 

0.13514643 

0.99880862 

0.13581328 

1.00799997 

0.320 

0.320 

0.13623553 

0.99073427 

0.13581449 

1 . OOOOOO^n 

0.000 

0.360 

0.11948419 

1.06435686 

0.12778250 

1.07199996 

0.040 

0.360 

0.12038714 

1.05629895 

0.12778273 

1.06399996 

0.080 

0.360 

0.12130390 

1.04824014 

0.12778300 

1 .05599996 

0.120 

0.360 

0.12223479 

1.04018044 

0.12778329 

1.04799997 

0.160 

0.360 

0.12318015 

1.03211979 

0.12778363 

1.03999997 

0.200 

0.360 

0.12414031 

1.02405819 

0.12778399 

1.03199997 

0.240 

0.360 

0.12511562 

1.01599560 

0.12778440 

1.02399997 

0.280 

0.360 

0.12610646 

1.00793202 

0.12778484 

1.01599997 
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0.320 

0.360 

0.12711318 

0.99986741 

0.12778531 

1.00799997 

0.360 

0.360 

0.12813620 

0.99180178 

0.12778584 

1 . ooooooon 

0.000 

0.400 

0.11112296 

1.07333873 

0*1 1976595 

1.07999996 

0.040 

0.400 

0.11195548 

1.06528872 

0.11976571 

1.07199996 

0.080 

0.400 

0.11280061 

1.05723796 

0.11976548 

1.06399996 

0.120 

0.400 

0.11365864 

1.04918645 

0.11976527 

1.05599997 

0.160  • 

-  0.400. 

0.11452988  ' 

1.04113411 

!o. 11976509 

1.04799997 

0.1200 

0.400 

0.11541464 

1.03308096 

0.11976492 

1.03999997 

0.240 

0.400 

0.11631322 

1 *02502697 

0.11976477 

1.03199997 

0.280 

0*400 

6.11722597 

1.01697215 

Qi 1 1 976465 

1 • 62399997 

0.320 

0.4oo 

6.11815321 

1.00891643 

0.11976455 

1.61599997 

0.360 

0.400 

0.11909531 

1.00085982 

0.11976448 

1.00799997 

0.400 

0.400 

0.12005261 

0.99280232 

0.11976444 

1.00000000 

0.000 

0.440 

0.10289846 

1.08224511 

0.11175607 

1.08799996 

0.040 

0.440 

0.10366283 

1.07420231 

0.11175545 

1.07999996 

0.080 

0.440 

0.10443868 

1.06615887 

0.11175484 

1.07199996 

0.120 

0.440 

0.10522626 

1.05811479 

0.11175424 

1.06399998 

0.160 

0.440 

0.10602586 

1.05007001 

0.11175365 

1.05599997 

0.200 

0.440 

0.10683774 

1.04202455 

0.11175306 

1.04799997 

0.240 

0.440 

0.10766219 

1.03397838 

0.11175249 

1.03999997 

0.280 

0.440 

0.10849951 

1.02593149 

0.11175193 

1.03199997 

0.320 

0.440 

0.10935000 

1.01788386 

0.11175138 

1.02399997 

0.360 

0.440 

0.11021399 

1.00983546 

0.11175084 

1.01599997 

0.400 

0.440 

0.11109178 

1.00178630 

0.11175032 

1.00799997 

0.440 

0.440 

0.11198372 

0.99373637 
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0.03461375 

1.15130992 

0.03986708 

1.15199998 

0.080 

0.800 

0.03485513 

1.14330512 

0.03986619 

1.14399997 

0.120 

0.800 

0.03509990 

1.13530026 

0.03986529 

1.13599997 

0.160 

0.800 

0.03534812 

1.12729533 

0.03986438 

1.12799999 

0.200 

0.800 

0.03559989 

1.11929033 

0.03986345 

1.11999998 
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240 

0.800 

280 

0.800 

320 

0.800 

360 

0.800 

400 

0.800 

440 

0.800 

480 

0.800 

520 

0.800 

560 

0.800 

600 

0.800 

640 

0.800 

680 

0.800 

720 

0.800 

760 

0.800 

800 

0.800 

000 

0.840 

040 

0.840 

080 

0.840 

120 

0.840 

160 

0.840 

200 

0.840 

240 

0.840 

280 

0.840 

320 

0.840 

360 

0.840 

400 

0.840 

440 

0.840 

480 

0.840 

520 

0.840 

560 

0.840 

600 

0.840 

640 

0.840 

680 

0.840 

720 

0.840 

760 

0.840 

800 

0.840 

840 

0.840 

000 

0.880 

040 

0.880 

080 

0.880 

120 

0.880 

160 

0.880 

200 

0.880 

240 

0.880 

280 

0.880 

320 

0.880 

360 

0.880 

400 

0.880 

440 

0.880 

480 

0.880 

520 

0.880 

560 

0.880 

600 

0.880 

640 

0.880 

680 

0.880 

720 

0.880 

0.03585527 

0.03611434 

0.03637718 

0.03664388 

0.03691452 

0.03718919 

0.03746798 

0.03775098 

0.03803829 

0.03833001 

0.03862625 

0.03892709 

0.03923267 

0.03954307 

0.03985844 

0.02730927 

0.02749704 

0.02768740 

0.02788043 

0.02807616 

0.02827467 

0.02847600 

0.02868022 

0.02888739 

0.02909757 

0.02931084 

0.02952726 

0.02974689 

0.02996983 

0.03019612 

0.03042587 

0.03065913 

0.03089600 

0.03113656 

0.03138090 

0.03162910 

0.03188127 

0.02034111 

0.02047998 

0.02062077 

0.02076351 

0.02090823 

0.02105499 

0.02120382 

0.02135477 

0.02150789 
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0.02198070 

0*02214296 

0.02230763 

0.02247477 

0.02264444 

0.02281668 

0.02299157 

0.02316916 


1.11128526 

1.10328011 

1.09527488 

1.08726958 

1.07926420 

1.07125875 

1.06325320 

1.05524760 

1.04724189 

1.03923609 
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1.01521818 

1.00721201 

0.99920575 
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1.03951121 
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1.07199999 
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1.05599999 
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1.03200000 


° '  .  r 

V  i  .  ?  ?  '  .  ■ 

'  '•  ,  r 

"  ^  A!  . 

W  [  /  i ’  •  ' 

.Off  .) 

3  rmaeo.o 

'?  V-V'.  „ 

'  ’«  : '  ..  I 

f.  •  o.r 

.  ■ 

r  .  r 

■ 

'  ■  I  , 

■  r  •  : 

0^.0 

)  v  . 

? 

r  ;  :  ' 

•  .-v.':  c,  . 

08.0 

£  '■  £0.0 

■  r  ■  ■  •  ;  *  f 

0  8*0 

■  ■  s  ■  ■  ■  • 

.  :<•  -  ■ 

- 

•  ' .  r 

■ 

T.-r*??,.  .  ,  i 

.  •  • 

:o  ••: 

T  -V  !  i  X  .  I 

■  •  .  ' 

Qth  ,.  ) 

.. 

'■■■■•.•  •'  8  .  . 

0^8.0 

;  .  ,  ps-  •  ,  - 

1  8  8.0 

■  0  6  •-  ,  1 

T  ?f9(  : s c * c 

!ri  0 . : 

919$? r  o , l 

■ 

r  - 

•r<  "0.0 

■  ■81  * 

'  '  '  o . 

M  *i 

V  >4-:  '  .. 

i'«  C?f'.04 

?  '  •  ’•>  . 

4  £ 

:  r  i.  o  • 

•  ■  .r 

C-  .  '  :  r 

...  .  , 

0  -v  —  ^ 

IV  •  .0 

'  •  **  •  • 

v£t38:  £■.  .C 

0-88.0 

I 

a 

..  ,  . .  . 

t^oseos  .0 

C  '8.0 

0  30*0 

■  ■ '  ,/ 

OSf  .0 

.  i 

*  '■  •  0  A ; : 

...  .  f 

’(  r  J 

o-^rvsr. 

*»  •• "  I  :•  ■  ' 

• r  ; ...  ■ 

~ r  •.  r.  i 

* 

,  .  •••  r  „ 

C  •  8  .  C 

SSC.C 

(  -  e .  0 

rr: '  o.r 

:  •  £!..’•  . 

(  vC 

■  ■  . 

•:r  '  0*  ■ 

•:  „ 

X  ■  £  - 

08  c.:) 

O^d.O 

•a  < 

109 


0.760 

0.880 

0.02334951 

1.02372086 

0.02390773 

1 • O24OO0OO 

0.800 

0.880 

0.02353270 

1.01571868 

0.02390701 

1 . 016OO0OO 

0.840 

0.880 

0.02371878 

1.00771646 

0.02390629 

1 . 00800000 

0.880 

0.880 

0.02390783 

0.99971421 

0.02390555 

1 • 000OO00O 

0.000 

0.920 

0.01346854 

1.18389257 

0.01594627 

1  .  1 839999Q 

0.040 

0.920 

0.01355986 

1.17589186 

0.01594591 

1 .17599998 

0.080 

0.920 

0.01365243 

1.16789111 

0.01594554 

1 .16799998 

0.120 

0.920 

0.01374627 

1 .15989038 

0.01594517 

1  .  1 5999998 

0.160 

0.920 

0.01384141 

1.15188963 

0.01594479 

1 .15199999 

0.200 

0.920 

0.01393787 

1.14388886 

0.01594441 

1 .14399999 

0.240 

0.920 

0.01403569 

1.13588808 

0.01594402 

1 .13599999 

0.280 

0.920 

0.01413489 

1.12788729 

0.01594362 

1.12799999 

0.320 

0.920 

0.01423550 

1.11988650 

0.01594322 

1  . 11999998 

0.360 

0.920 

0.01433756 

1.11188568 

0.01594282 

1.11199998 

0.400 

0.920 

0.01444109 

1 .10388486 

0.01594241 

1 .10399998 

0.440 

0.920 

0.01454612 

1.09588402 

0.01594199 

1 . 09599998 

0.480 

0.920 

0.01465270 

1.08788317 

0.01594156 

1.08799998 

0.520 

0.920 

0.01476085 

1.07988233 

0.01594113 

1 .07999999 

0.560 

0.920 

0.01487060 

1.07188145 

0.01594070 

1 • 07199909 

0.600 

0.920 

0.01498200 

1.06388057 

0 . o l  5  04°2  5 

1 . 06399999 

0  •  640 

0.920 

0.01509509 

1.05587966 

0.01593980 

1 .05509999 

0.680 

0.920 

0.01520989 

1.04787875 

0.01593935 

1 . 04799999 

0.720 

0.920 

0.01532645 

1.03987784 

0.01593889 

1 . ^4000000 

0.760 

0.920 

0.01544481 

1.03187690 

0.01593841 

1 . 03200000 

0.800 

0.920 

0.01556502 

1.02387595 

0.01593794 

1 . 0240O0OO 

0.840 

0.920 

0.01568711 

1.01587498 

0.01593745 

1 .01600000 

0.880 

0.920 

0.01581113 

1.00787399 

0.01593696 

1 . 008000OO 

0.920 

0.920 

0.01593713 

0.99987300 

0.01593646 

1 • OOOOooOO 

0.000 

0.960 

0.00668900 

1.19197331 

0.00797323 

1 .19199998 

0.040 

0.960 

0.00673404 

1.18397311 

0.00797305 

1 .18399998 

0.080 

0.960 

0.00677970 

1.17597295 

0.00797287 

1 . 17599998 

0.120 

0.960 

0.00682598 

1.16797275 

0.00797268 

1.16799908 

0.160 

0.960 

0.00687289 

1.15997259 

0.00797249 

1.15999999 

0.200 

0.960 

0.00692046 

1.15197239 

0.00797230 

1.15109999 

0.240 

0.960 

0.00696868 

1.14397220 

0.00797211 

1 .14399999 

0.280 

0.960 

0.00701759 

1.13597200 

0.00797191 

1 .13599999 

0.320 

0.960 

0.00706718 

1.12797181 

0.00797171 

1.12799999 

0.360 

0.960 

0.00711748 

1.11997161 

0.00797151 

1.11999998 

0.400 

0.960 

0.00716850 

1.11197140 

0.00797131 

1.11199998 

0.440 

0.960 

0.00722026 

1.10397119 
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1 .10399998 

0.480 

0.960 
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1.09597098 
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1.09599998 

0.520 

0.960 

0.00732606 

1.08797079 
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1 .08799999 

0.560 

0.960 
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1.07997056 

0.00797046 

1 . 07909999 

0.600 

0.960 

0.00743500 

1.07197035 

0.00797024 

1 • 07199999 

0.640 

0.960 

0.00749069 

1.06397013 

0 .00797002 

1 .06309999 

0.680 

0.960 

0.00754722 

1.05596990 

0.00796979 

1 .05599999 

0.720 

0.960 

0.00760462 

1.04796969 

0.00796956 

1 . 048OO0O0 

0.760 

0.960 

0.00766289 

1.03996945 

0.00796933 

1 . 04000000 

0.800 

0.960 

0.00772206 

1.03196921 

0.00796909 

1 . 03200000 

0.840 

0.960 

0.00778216 

1.02396897 

0.00796885 

1 .02400000 

0.880 

0.960 

0.00784320 

1.01596873 

0.00796861 

1.01600000 

0.920 

0.960 

0.00790520 

1.00796851 

0.00796836 

1 • 00800001 

0.960 

0.960 

0.00796819 

0.99996825 

0.00796811 

1 • 0000^000 

-0.000 

1.000 

0. 09000000 

1 .19999998 
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0.040 
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0.120 

1.000 

0.00000000 

1 . 1 7599998 

0 .  9*90090* 

1 .17599998 

0.160 

1.000 

0. 0^00^000 

1 .16799999 

0 ,000^00°^ 

1 . If 700000 

0*200 

1 . 000 

0.00000000 

1.15999999 

0  .  ****09*0 

1 . 1 5999909 

0.240 

1.000 

0. OOOOHOOO 

1.15199999 

0  .9**00900 

1  •  1 51QQ990 

0.280 

1.000 

0 .00000000 

1.14399999 

0.909**0*9 

1 .14399999 

0.320 

1.000 

0.00000000 

1 .13599999 

0  .oooooooo 

1 . 13599999 

0.360 

1.000 

0.00000000 

1.12799999 

0 .OOOOOOOO 

1 .12799999 

0.400 

1.000 

0.00000000 

1.11999998 

0. OOOOOOOO 

1.11999998 

0.440 

1.000 

0.00000000 

1.11199998 

0.00000000 

1.11199998 

0.480 

1.000 

o. 'oooooooo 

1.10399998 

0 • oooooooo 

1 .10399998 

0.520 

1.000 

0 • 000000^0 

1 .09600000 

* . *00*0000 

1 • 

0.560 

1  •  poo 

0. OOOOOOOO 

l .0879999Q 

0. OOOOOOOO 

1 • 0R79Q999 

0.600 

1.000 

0.00000000 

} .07099990 

0 . 0  00090*0 

1 . *7000009 
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i. 67199999 
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0.680 

1 . 000 
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0  •  0  9000900 

1 , 06399999 

0.720 

1 . 000 

0. 00000000 

1 .05600000 

0 .  0*9*00*9 

1 . o56n Aoo* 

0.760 

1. 000 

0  •  oooooooo 

1 . *480*00* 
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1 . *48000*0 

0.800 

1. 000 

0.00000000 

1 . 04000O0O 

n  • 

1 . 94*00909 

0.840 

1.000 

0.00000000 

1 .03200000 
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0. oooooooo 
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1 . *240**** 
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1.000 
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1 .01600001 
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COMPUTER  PROGRAM  P 


READ  (5*100)  E  *  H  »  N 

HH=H*H 

G=  0 • 96 

WRITE  (6*101) 

CC  =  (  1 . -G ) *16. /3. 141 5 9265/3 ,141 5926 5 
CB=2.*Ct/3. 14159265 
TQ  =  0  * 

J=  N 

Y=l.+H/2* 

10  JsJ-1 
Y=  V-H 
X--H/2. 

DO  30  1  =  1, J 
X=  X  +  H 

IE  ( I.LT. J)  GO  TO  12 

YY=Y 

Y= Y+H/6 • 

X=  X-H/6 • 

HH=0 • 5*HH 
12  SGa-1.0 

ww=o , 

WX  =  0. 

WY  =  0  • 

DO  20  N 1  =  1 ,8 
FN=N I “1 

GG=(EN+0.5)*3.14159265 

ANC=0.5*(EXP(GG)+EXP(-GG) )*( l.+2.*FN)*(l .+2.*FN) 
ANB  =  ANC*  (  1  •  4-2 •  *FN  ) 

GGX±GG*X 

GGY=GG*Y 

15  EHY=0.5*( EXP(GGY)-EXP(-GGY) ) 

CH Y=EXP ( -GGY ) +SHY 
SG=-1 • *SG 

WW=WW+SG*SHY*S  I  N  (  GGX  )  /A.MB 
WX*WX+$G*SHY*COS  (  GGX  )  /ANC 
20  WY=WY+SG*CHY*SIN(GGX)/ANC 
EX=-2.*G*( Y-X)-CC*WX 
EY=2.*X*( l.-G)-CC*WY 
Wte-G*X*(Y“X)*f  ( l.-G)*X*Y~C»*WW 
S=EX*EX+EY*EY 
SS=SQRT(S) 

B=E+1./SS 

$X=B*EX 

SY=B*EY 

$T=SQRT(SX*SX+SY*SY) 

TQ=TO+ ( SY*X-SX*Y ) *HH 
30  WRITE  (6*102)  X * Y , W , SX , SY , ST 
WRITE  (6*103) 

32  Y=  YY 

HH=2.*HH 

IF  (J.GT.l)  GO  TO  10 
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TQQ=TQ*8. 


WRITE 

WRITE 

STOP 

(6*104)  TQQ 
(6*105)  TQ 

100 

FORMAL 

(F10.3,F9. 

3*16) 

101 

FORMAT 

1  sv 

( 74H 

X  Y 

ST/  ) 

102 
1  08 

FORMAT 

FORMAT 

(8X*2F8.3* 

(1HJ) 

2F14.8 ,2F13.3 ) 

104 

FORMAT 
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X 

y 

w 

s 

xz 

s 

yz 

T 

0.030 

0.990 

-0.02805391 

-1.18963275 

0.00005464 

1 .18963275 

0.070 

0.990 

-0.06277041 

-1.18194842 

0.00009264 

1 .18194842 

0.110 

0.990 

-0.09441099 

-1.17425154 

0.00007282 

1 .17425154 

0.150 

0.990 

-0.12296902 

-1.16653510 

0.00002812 

1 .16653510 

0.190 

0.990 

-0.14843657 

-1.15879999 

0.00002084 

1 .15879999 

0.230 

0.990 

-0.17080761 

-1.15105420 

0.00008623 

1 .15105420 

0.270 

0.990 

-0.19007957 

-1.14330528 

0.00019541 

1 .14330528 

0.310 

0.990 

-0.20625146 

-1.13555298 

0.00027515 

1 .13555301 

0.350 

0.990 

-0.21932039 

-1.12778860 

0.00027356 

1.12778863 

0.390 

0.990 

-0.22927941 

-1 .12000223 

0.00021604 

1 .12000225 

0.430 

0.990 

-0.23611893 

-1.11219161 

0.00019626 

1.11219163 

0.470 

0.990 

-0.23983050 

-1.10436429 

0.00029875 

1 .10436432 

0.510 

0.990 

-0.24040958 

-1.09653037 

0.00051326 

1.09653849 

0.550 

0.990 

-0.23785440 

-1.08869212 

0.00072228 

1.08869236 

0.590 

0.990 

-0.23216159 

-1.08083954 

0.0007937O 

1.08083982 

0.630 

0.990 

-0.22332221 

-1.07295688 

0.00072014 

1.07295711 

0.670 

0.990 

-0.21132207 

-1.06503556 

0.00068638 

1.06503577 

0.710 

0.990 

-0.19614656 

-1.05708213 

0.00097537 

1.05708257 

0.750 

0.990 

-0.17778601 

-1.04911182 

0.00173759 

1.04911325 

0.790 

0.990 

-0.15623607 

-1.04113145 

0.00277893 

1.04113516 

0.830 

0.990 

-0.13149073 

-1.03312508 

0.00358937 

1.03313131 

0.870 

0.990 

-0.10353220 

-1.02505423 

0.00382050 

1.02506133 

0.910 

0.990 

-0.07232507 

-1.01687207 

0.00468588 

1.01688264 

0.950 

0.990 

-0.03782147 

-1.00846590 

0.01523320 

1.00858093 

0.987 

0.993 

-0.00653873 

-0.98996323 

0.15190378 

1.00154976 

0.030 

0.950 

-0.02685800 

-1.18164528 

0.00008094 

1 .18164528 

0.070 

0.950 

-0.05997844 

-1.17395524 

0.00018129 

1.17395524 

0.110 

0.950 

-0.09002070 

-1.16625380 

0.00026745 

1 .16625382 

0.150 

0.950 

-0.11697963 

-1.15853834 

0.00035366 

1.15853838 

0.190 

0.950 

-0.14084956 

-1.15080904 

0.00046569 

1.15080913 

0.230 

0.950 

-0.16162549 

-1.14306874 

0.00061921 

1 .14306890 

0.270 

0.950 

-0.17930357 

-1.13531991 

0.00080583 

1 .13532828 

0.310 

0.950 

-0.19388046 

-1.12756213 

0.00100123 

1 .12756258 

0.350 

0.950 

-0.20535195 

-1.11979183 

0.00119090 

1 .11979246 

0.39  0 

0.950 

-0.21371220 

-1.11200509 

0.00139206 

1.112 0° 596 

0.430 

0.950 

-0.21895431 

-1.10420073 

0.00165007 

1.10420195 

0.470 

0.950 

-0.22107167 

-1.09638099 

0.00200862 

1.09638282 

0.510 

0.950 

-0.22005887 

-1.08854882 

0.00247802 

1.08855164 

0.550 

0.950 

-0.21591121 

-1.08070397 

0*00303395 

1.08070822 

0.590 

0.950 

-0.20862295 

-1.07284153 

0.00365969 

1.07284776 

0.630 

0.950 

-0.19818592 

-1.06495444 

0.00440892 

1.06496356 

0.670 

0.950 

-0.18458968 

-1.05703781 

0.00544428 

1.05705183 

0.710 

0.950 

-0.16782330 

-1 . 04909057 

0. 007O2289 

1.04911407 

0.750 

0.950 

-0. 14787714 

-1.04111129 

0.00945957 

1.04115426 

0.790 

0.950 

-0.12474235 

-1.03308749 

0.01319230 

1.03317171 

0.830 

0.950 

-0.09840790 

-1.02497515 

0 • 01 92  c  690 

1.02515602 

0.870 

0.950 

-0.06885661 

-1.01660764 

n  •0313r>6r' 3 

1 .01 78RQ54 

0.910 

0.950 

-0.03606374 

-1.00665608 

0.06817853 

1.08896221 

0.947 

0.953 

-0.00626325 

-0.95153497 

0.31455589 

1.00217972 

0.030 

0.910 

-0.02566421 

-1.17366558 

0.00012729 

1 .17366558 

0.070 

0.910 

-0.05719239 

-1.16597355 

0.00030372 

1  .  16597358 

0.110 

0.910 

-0.08564164 

-1 .15827062 

0.00048792 

1.15827071 

0.150 

0.910 

-0.11100738 

-1.15055574 

0.00068819 

1 .15055593 
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0.190 

0.910 

-0.13328481 

-1.14282889 

0.00091785 

1.14282925 

0.230 

0.910 

-C. 15246930 

-1.13509090 

0.00118715 

1 .13509151 

0.270 

0.910 

-0.16855661 

-1.12734246 

0.00149827 

1.12734344 

0.310 

0.910 

-0.18154256 

-1.11958307 

0.00184891 

1 .11958459 

0.350 

0.910 

-0.19142253 

-1.11181111 

0.00224287 

1.11181336 

0.390 

0.910 

-0.19819119 

-1  •  10402470 

n .  00269932 

1 . 1 o4o2  s  on 

0*430 

0.910 

-0.20184270 

-1.09622292 

0.00325243 

1.09622773 

0.470 

0.910 

-0.20237114 

-1.08840593 

0.00394171 

1.08841305 

0*510 

0.910 

-0.19977086 

-1.08057382 

0.00480295 

1 .08058449 

0.550 

0.910 

-0.19403616 

-1.07272513 

0.00587385 

1.07274121 

0.590 

0.910 

-0.18516067 

-1.06485620 

0.00722106 

1.06488068 

0.630 

0.910 

-0.17313677 

-1.05696146 

0.00898366 

1.05699963 

0.670 

0.910 

-0.15795569 

-1.04903373 

0.01142522 

1.04909594 

0.710 

0.910 

-0.13960808 

-1.04106107 

0.01500901 

1.04116925 

0.750 

0.910 

-0.11808454 

-1.03301465 

0.02058808 

1.03321978 

0.79  0 

0.910 

-0.09337520 

-1.02480581 

0.03004452 

1.02524612 

0.830 

0.910 

-0.06546848 

-1.01607021 

0.04887267 

] .01724490 

0.870 

0.910 

-0.03434979 

-1.00402070 

0.10232536 

1 .000221 49 

0.907 

0.913 

-0.00597513 

-0.92261163 

0.39246111 

1.00261555 

0.030 

0.870 

-0.02447318 

-1.16569536 

0.00018377 

1 .16569536 

0.070 

0.870 

-0.05441318 

-1.15800276 

0.00044554 

1 .15800283 

0.110 

0.870 

-0.08127405 

-1 .15029950 

0.00072981 

1 .15029972 

0.150 

0.870 

-0.10505149 

-1.14258517 

0.00104412 

1.14258565 

0.190 

0.870 

-0.12574105 

-1.13485959 

0.00139896 

1 .13486045 

0.230 

0.870 

-0.14333832 

-1.12712286 

0.00180501 

1.12712430 

0.270 

0.870 

-0.15783894 

-1.11937503 

0.00227178 

1.11937732 

0.310 

0.870 

-0.16923850 

“*1.11161552 

0.00280968 

1.11161905 

0.350 

0.870 

-0.17753231 

-1.10384337 

0.00343519 

1.10384871 

0.390 

0.870 

-0.18271535 

-1.09605738 

0.00417598 

1.09606532 

0.430 

0.870 

-0.18478221 

-1.08825644 

0.00507316 

1.08826826 

0.470 

0.870 

-0.18372735 

-1.08043965 

0. 00618146 

1.08045732 

0.510 

0.870 

-0.17954514 

-1.07260550 

0.00757235 

1.07263223 

0.550 

0.870 

-0.17222968 

-1.06475104 

0.00934823 

1.06479208 

0.590 

0.870 

-0.16177458 

“1 .05687089 

0.01167516 

1.05693536 

0.630 

0.870 

-0.14817272 

-1.04895535 

0.01484419 

1.04906037 

0.670 

0.870 

-0.13141624 

-1.04098519 

0.01938896 

1.04116573 

0.710 

0.870 

-0.11149676 

-1.03291487 

0.02635706 

1.03325110 

0.750 

0.870 

-0.08840546 

-1.02460853 

0.03810139 

1.02531670 

0.790 

0.870 

-0.06213288 

-1.01551657 

0.06129912 

1.01736497 

0.830 

0.870 

-0.03266849 

-1.00160877 

0.12527246 

1.00941236 

0.867 

0.873 

-0.00569148 

-0.90373711 

0.43488313 

] .00292^7* 

0.030 

0.830 

-0.02328513 

-1.15773518 

0.00024682 

1.15773519 

0.07  0 

0.830 

-0.05164104 

-1 .15004250 

0.0006^194 

1 .15^04265 

0.110 

0.830 

-0.07691788 

-1 .14233939 

0.00099277 

1 .14233982 

0.150 

0.830 

-0.09911148 

-1.13462562 

0.00142844 

1 .13462652 

0.190 

0.830 

-0.11821761 

-1.12690094 

0.00192079 

1 .12690258 

0.230 

0.830 

-0.13423198 

-1.11916518 

0.00248386 

1.11916793 

0.270 

0.830 

-0.14715028 

-1.11141803 

0.00313421 

1 .11142244 

0.310 

0.830 

-0.15696810 

-1.10365892 

0.00389289 

1 .10366578 

0.350 

0.830 

-0.16368084 

-1.09588691 

0.00478920 

1.09589737 

0.390 

0.830 

-0.16728364 

-1.08810085 

0.00586544 

1.08811664 

0.430 

0.8  30 

-0.16777142 

-1.08029919 

0.00718226 

1.08032307 

0.470 

0.830 

-0.16513889 

-1.07247990 

0.00882605 

1.07251621 

0.510 

0.830 

-0.15938058 

-1 .06463958 

0.91^92342 

1 .0646Q561 

0.550 

0.830 

-0.15049070 

-1.05677215 

0.01367107 

1 .05686058 

0.590 

0.830 

-0.13846312 

-1.04886591 

0.01739763 

1.04901018 

0.630 

0.830 

-0.12329118 

-1.04089620 

0.02269515 

1.04114358 

0.670 

0.830 

-0.10496773 

-1.03280331 

0.03073202 

1.03326043 
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0.710 

0.830 

-0.08348517 

-1.02441082 

0.04414581 

1.02536158 

0.750 

0.830 

-0.05883547 

-1.01501539 

0 . 07036674 

1.01745158 

0.79  0 

0.830 

-0.03101011 

-0.99961906 

0.14128029 

1.00955355 

0.827 

0.833 

-0.00541296 

-0.89151434 

0.45991788 

1.00315615 

0.030 

0.790 

-0.02210008 

-1.14978502 

0 . 00031539 

1 .14978505 

0.070 

0.790 

-0.04887597 

-1.14209249 

^.00^77194 

1 .14209274 

0.110 

0.790 

-0.07257293 

-1.13438974 

0.00127834 

1 . 13439H45 

0.150 

0.790 

-0.09318694 

-1.12667655 

0.00184670 

1  .  12667806 

0.190 

0.790 

-0.11071391 

-1.11895267 

0.00249248 

1  . 1 1895543 

0.230 

0.790 

-0.12514969 

-1.11121773 

0.00323520 

1.11122244 

0.270 

0.790 

-0.13649008 

-1.10347123 

0.00409980 

1 .10347883 

0.310 

0.790 

-0.14473076 

-1.09571242 

0.00511937 

1.09572437 

0.350 

0.790 

-0.14986725 

-1.08794018 

0.00633954 

1.08795864 

0.390 

0.790 

-0.15189492 

-1.08015286 

0.00782512 

1.08018121 

0.430 

0.790 

-0.15080891 

-1.07234805 

0.00967016 

1.07239165 

0.470 

0.790 

-0.14660419 

-1.06452177 

0.01201490 

1.06458957 

0.510 

0.790 

-0.13927551 

-1.05666701 

0.01507750 

1.05677457 

0.550 

0.790 

-0.12881736 

-1.04877015 

0.01921797 

1.04894621 

0.590 

0.790 

-0.11522394 

-1.04080215 

0.02507687 

1.04110419 

0.630 

0.790 

-0 . 09848909 

-1.03269198 

0.03391057 

1 .03324858 

0,670 

0.790 

-0.07860628 

-1.02423089 

^.04854749 

1 . 02  538^8^ 

0.710 

0.790 

-0.05556865 

-1.01459540 

0.0769-1  380 

1*01750653 

0.750 

0.790 

-0.02936903 

-0.99806504 

0.  ]  5257^42 

1.00965309 

0.787 

0.793 

-0.00513803 

-0.88347022 

0.47553545 

1.00332127 

0.03  0 

0.750 

-0.02091798 

-1.14184471 

0.00038944 

1 .14184477 

0.070 

0.750 

-0.04611782 

-1.13415244 

0.00095625 

1.13415284 

0.110 

0.750 

-0.06823893 

-1.12645013 

0.00158927 

1 .12645125 

0.150 

0.750 

-0.08727743 

-1.11873757 

0.00230479 

1 .11873994 

0.190 

0.750 

-0.10322938 

-1.11101439 

0.00312372 

1.11101877 

0.230 

0.750 

-o.  11609079 

-1.10328011 

0  •  00407328 

1.10328762 

0.270 

0.750 

-0.12585757 

-1.09553398 

0.00518946 

1.09554627 

0.310 

0.750 

-0.13252556 

-1 . 08777494 

0.00652124 

1 .08779448 

0.350 

0.750 

-0.13609046 

-1 .08000132 

0.00813729 

1 .08^03  1 06 

0.390 

0.750 

-0.13654784 

-1.07221052 

0.01013690 

1.07225842 

0.430 

0.750 

-0.13389308 

-1.06439813 

0.01266859 

1.06447352 

0.470 

0.750 

-0.12812140 

-1.05655643 

0.01596397 

1.05667702 

0.510 

0.750 

-0.11922785 

-1.04867021 

0.02040457 

1.04886870 

0.550 

0.750 

-0.10720723 

-1.04070703 

0.02666610 

1.04104860 

0.590 

0.750 

-0.09205414 

-1.03258749 

0.03606773 

1.03321720 

0.630 

0.750 

-0.07376290 

-1.02407886 

0.05156840 

1.02537642 

0.670 

0.750 

-0.05232759 

-1.01426979 

0.08142275 

1 .01 753270 

0.710 

0.750 

-0.02774206 

-0.99692785 

0.16019258 

1 .00071619 

0.747 

0.753 

-0.00486574 

-0.87820502 

0.48542501 

1 .00343485 

0.030 

0.710 

-0.01973875 

-1.13391399 

0.00046941 

1.13391408 

0.07  0 

0.710 

-0.04336638 

-1 .12622198 

0.00115653 

1.12622256 

0.110 

0.710 

-0.06391554 

-1.11852015 

0.00192944 

1.11852180 

0.150 

0.710 

-0.08138247 

-1 .11080817 

0.00280998 

1.11081172 

0.190 

0.710 

-0.09576339 

-1  .  10308559 

0.00382661 

1.10309222 

0.230 

0.710 

-0.10705448 

-1.09535166 

0.00501711 

1 • 095363  1  c 

0.270 

0.710 

-0.11525182 

-1.08760534 

0.00643280 

1.08762437 

0.310 

0.710 

-0.12035142 

-1.07984496 

0.00814525 

1.07987567 

0.350 

0.710 

-0.12234917 

-1.07206779 

0.01025751 

1.07211685 

0.390 

0.710 

-0.12124086 

-1.06426930 

0.01292339 

1.06434775 

0.430 

0.710 

-0.11702213 

-1.05644117 

0.01638268 

1.05656819 

^•470 

0.710 

-0.10968846 

-1.04856725 

0.02103013 

1 .04877PH 

0.510 

0.710 

-0.09923521 

-1.04061274 

0.02756376 

1.04097773 

0.550 

0.710 

-0.08565754 

-1.03249271 

0.03734312 

1.03316779 

0.590 

0.710 

-0.06895044 

-1 .02395859 

0.05341008 

1.02535058 
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0.630 

0.710 

-0.04910874 

-1.01404149 

0.08422368 

1.01753318 

0.670 

0.710 

-0.02612708 

-0.99617953 

0.16497271 

1.00974732 

0.707 

0.713 

-0.00459561 

-0.87495660 

0.49139864 

1 .00350469 

0.030 

0.670 

-0.01856229 

-1.12599249 

0.00055609 

1.12599262 

0.070 

0.670 

-0.04062143 

-1.11830078 

0.00137525 

1 .11830161 

0.110 

0.670 

-0.05960235 

-1.11059939 

0.00230409 

1 .11060178 

0.150 

0.670 

-0.07550149 

-1.10288792 

0.00337183 

1 .10289306 

0.190 

0.670 

-0.08831523 

-1.09516564 

0.00461720 

1.09517537 

0.230 

0.670 

-0.09803991 

-1.08743156 

0.00609284 

1.08744862 

0.270 

0.670 

-0.10467180 

-1 .07968400 

0.00787201 

1.07971268 

0.310 

0.670 

-0.1082071 1 

-1.07192024 

0. 01 005998 

1.071 96745 

0.350 

0.670 

-0.10864193 

-1.06413565 

0.01281355 

1.06421279 

0.390 

0.670 

-0.10597229 

-1.05632174 

0.01637686 

1.05644866 

0.430 

0.670 

-0.10019408 

-1.04846182 

0.02115151 

1.04867514 

0.470 

0.670 

-0.09130309 

-1.04051999 

0.02784703 

1.04089255 

0.510 

0.670 

-0.07929496 

-1.03240841 

0.03784387 

1.03310176 

0.550 

0.670 

-0.06416524 

-1.02387031 

0 . 05422660 

1.02530530 

0.590 

0.670 

-0.04590933 

-1.01390704 

0.08555839 

1.01751054 

0.630 

0.670 

-0.02452251 

-0.99578736 

0.16734055 

1.00975012 

0.667 

0.673 

-0.00432732 

-0.87328900 

0.49442115 

1.00353671 

0.030 

0.630 

-0.01738849 

-1 . 1 1807985 

0.00065057 

1 • 1 1 808O07 

0.070 

0.630 

-0.03788267 

-1.11038841 

0.00161579 

1.11 038957 

0 . 1 1  0 

0.6  30 

-0.05529895 

-1.10268740 

0.00272028 

1 .10269073 

0.150 

0.630 

-0.06963301 

-1 .09497620 

0.00400320 

1.09498351 

0.190 

0.630 

-0.08088413 

-1.08725382 

0.00551762 

1.08726781 

0.230 

0.630 

-0.08904612 

-1.07951864 

0.00733741 

1.07954357 

0.270 

0.630 

-0.09411636 

-1.07176807 

0.00956857 

1 .07181078 

0.310 

0.630 

-0.09609125 

-1.06399744 

0.01236880 

1.06406933 

0.350 

0.630 

-0.09496714 

-1.05619832 

0.01598339 

1.05631925 

0.390 

0.630 

-0.09074027 

-1.04835411 

0.02081558 

1.04856073 

0.430 

0.630 

-0.08340680 

-1 .04042872 

0.02757732 

1.04079413 

o  .470 

0.630 

-0.07296282 

-1.0323340 2 

0.03765329 

1.03302046 

0.510 

0.630 

-0.05940429 

-1.02381202 

0.05413610 

1.02524229 

0.550 

0.630 

-0.04272712 

-1.01385936 

0.08560752 

1 .01746717 

0.59  0 

0.630 

-0.02292710 

-0.99571710 

0.16762368 

1.00972780 

0.627 

0.633 

-0.00406066 

-0.87293921 

0.49503621 

1 •00353561 

0.030 

0.590 

-0.01621721 

-1.11017562 

0.00075433 

l.ll0l7587 

0.070 

0.590 

-0.03514983 

-1.10248440 

0.00188271 

1 .10248599 

0.110 

0.590 

-0.05100485 

-1.09478363 

0.00318761 

1.09478826 

0.150 

0. 590 

-0.06377907 

-1.08707238 

0.00472189 

1.08708263 

0.190 

0.590 

-0.07346923 

-1.07934913 

0 • 00655906 

1.07936905 

0.230 

0.590 

-0.08007207 

-1.07161138 

0.00880452 

1.07164755 

0.270 

0.590 

-0.08358425 

-1.06385472 

0.01161498 

1.06391811 

0.310 

0.590 

-0 . 08400240 

-1.05607095 

0.01523403 

1.05618081 

0.350 

0.590 

-0.08132309 

-1.04824389 

0.02006193 

1 .04843583 

0  i  390 

0.590 

-3.07554282 

-1 . ^4033838 

O.n268o528 

1.04063763 

0.430 

0.590 

-0.06665803 

-1.03226816 

0.03683829 

1 .07202^2- 

0.470 

0.590 

-0.05466507 

-1.02378033 

0.05323100 

1.02516324 

0.510 

0.590 

-0.03956026 

-1.01388932 

0.08450826 

1.01740513 

0.550 

0.590 

-0.02133983 

-0.99593463 

0.16605439 

1 .00968304 

0.587 

0.593 

-0.00379545 

-0.87374780 

0.49354573 

1. 00350515 

o.o30 

0*550 

-0.01504833 

-1.10227935 

0 • 00086936 

1.1o227969 

0.070 

0*550 

-0.03242258 

-1.09458829 

0.00218222 

1.09459045 

0.110 

0.550 

-0.04671957 

-1 . 08688750 

0.00371933 

1.08689387 

0.150 

0.550 

-0.05793628 

-1.07917562 

0.00555300 

1.07918989 

0.190 

0.550 

-0.36606967 

-1.07145033 

o .00778686 

1.07147862 

0.230 

0.550 

-0.07111666 

-1 .06370748 

0.01057490 

1 .06376004 

^.270 

0.550 

-0.0730741 5 

-1 .05593934 

o. 01415644 

1.05603422 
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0.310 

0.550 

-0.07193899 

-1.04813064 

0.01892465 

1.04930147 

0.350 

0.550 

-0.06770798 

-1.04024790 

0.02 557370 

1.04056221 

0.390 

0.550 

-0.06037789 

-1.03220875 

0.03545415 

1.03281745 

0.43  0 

0.550 

-0.04994540 

-1.02377092 

0.05158567 

1.02506974 

0.470 

0.550 

-0.03640718 

-1.01393646 

0.08236700 

1.01732632 

0.510 

0.550 

-0.01975984 

-0.99640640 

0.16279829 

1.00961823 

0.547 

0.553 

-0.00353152 

-0.87562488 

0.49009199 

1 .003448  5  9 

0.030 

0.510 

-0.01388170 

-1.09439055 

0.0O099836 

1.09439100 

0.070 

0.510 

-0.02970058 

-1.08669951 

0.00252289 

1.08670243 

0.110 

0.510 

-0.04244255 

-1.07899837 

0.00433418 

1.07900707 

0.150 

0.510 

-0.05210480 

-1.07128499 

0.00653302 

1.07139491 

0.190 

0.510 

-0.05868449 

-1.06355564 

0.00926916 

1.06359602 

0.230 

0.510 

-0.06217875 

-1.05580326 

0.01277527 

1 . 0*  588054 

0.270 

0.510 

-0.06258470 

-1.04801371 

0.01743366 

1 • 0  4  ° 1 *070 

0.310 

0.510 

-0.05989943 

-1.04015599 

0.02391958 

1.04043097 

0.350 

0.510 

-0.05411997 

-1.03215332 

0.03354777 

1.03269836 

0.390 

0.510 

-0.04524333 

-1.02377875 

0.04926170 

1.02496324 

0.43  0 

0.510 

-0.03326650 

-1.01413926 

0.07926841 

1.01723247 

0.470 

0.510 

-0.01818640 

-0.99709898 

0.15797365 

1.00953555 

0.507 

0.513 

-0.00326874 

-0.87853109 

0.48469768 

1.00336869 

0.030 

0.470 

-0.01271716 

-1.08650875 

0.00114504 

1.08650934 

0.070 

0.470 

-0.02698349 

-1.07881753 

0.00291680 

1.07982146 

0.110 

0.470 

-0.03817326 

-1.07111540 

0.00505936 

1.07112734 

0.150 

0.470 

-0  .  ''>4628388 

-1 . 06339907 

o. 00771679 

1.06342707 

0.190 

0.470 

-0.05131272 

-1.05566224 

0.01111295 

1.05572073 

0.230 

0.470 

-C. 05325714 

-1.04789225 

0.01561581 

1.04900859 

0.270 

0.470 

-0.05211448 

-1 .04006110 

0.02187566 

1.04029113 

0.310 

0.470 

-0.04788204 

-1 .03209897 

0.03115997 

1.03256922 

0.350 

0.470 

-0.04055712 

-1.02379825 

0.04631142 

1.02484515 

0.390 

0.470 

-0.03013697 

-1.01433548 

0.07528167 

1.01712526 

0.430 

0.470 

-0.01661885 

-0.99797838 

0.15166466 

1.00943697 

0.467 

0.473 

-0.00300699 

-0.88246611 

0.47728410 

1.00326792 

0.030 

0.430 

-0.01155458 

-1.07863338 

0.00131465 

1.07863417 

0.070 

0.430 

-0.02427094 

-1.07094168 

0.00338152 

1.07094702 

0.110 

0.430 

-0.03391113 

-1.06323763 

0.00593577 

1.06325420 

0.150 

0.430 

-0.04047274 

-1.05551590 

0.00919043 

1.05555590 

0.190 

0.430 

-0.04395338 

-1.04776537 

0.01349585 

1.04785227 

0.230 

0.430 

-0.04435061 

-1.03996156 

0.01947167 

1.04014383 

0.270 

0.430 

-0.04166202 

-1.03204266 

0.02832723 

1.03243134 

0.310 

0.430 

-0.03588513 

-1.02382348 

0.04278065 

1.02471688 

0.350 

0.430 

-0.02701748 

-1.01456214 

0.07046502 

1.01700621 

°  •  39  0 

0.430 

-0.01505660 

-0.99900934 

0.14393087 

1.00932440 

0.427 

0.433 

-0.00274618 

-0.88746141 

0.46767433 

1.00314855 

0.030 

0.390 

-0.01039378 

-1.07076390 

0.00151482 

1.07076495 

0.070 

0.390 

-0.02156257 

-1.06307123 

0.00394347 

1.06307853 

0.110 

0.390 

-0.02965557 

-1.05536372 

0.00702769 

1.05538712 

0.150 

0.390 

-0.03467058 

-1.04763212 

0.011 09709 

1 .047^9089 

0.190 

0.390 

-0.03660543 

-1.03985558 

0.01673530 

] .03999024 

0.230 

0.390 

-0.03545792 

-1*03198117 

0.02508305 

1.03228594 

0*270 

0.390 

-0.03122584 

-1.02384822 

0.03871073 

1.0245797^ 

0.310 

0.390 

-0.02390698 

-1.01480565 

0.06486918 

1.01687684 

0.350 

0.390 

-0.01349910 

-1.00015446 

0.13481420 

1.00919958 

0.387 

0.393 

-0.00248621 

-0.89357514 

0.45558540 

1.00301274 

C.030 

0.350 

-0.00923461 

-1.06289973 

0.00175700 

1.06290118 

0.070 

0.350 

-0.01885801 

-1.05520531 

0.00464419 

1.05521552 

0.110 

0.350 

-0.02540598 

-1.04749153 

0.00844204 

1.04752553 
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0.150 

0.350 

-0.02887658 

-1.03974141 

0.01369298 

1.03983156 

0.190 

0.350 

“0.02926784 

-1.03191114 

0.02145908 

1.03213424 

0.230 

0.350 

“0.02657779 

-1 .0238661 0 

0.03414014 

1.02443511 

0.270 

0.350 

“0.02080445 

“1 .01505197 

0.05854004 

1.01673862 

0.310 

0.350 

-0.01194584 

-1.00137364 

0.12434439 

1 .00^06 42^ 

0.347 

0.353 

-0.00222699 

-0.90088760 

0.44060739 

1.00286257 

0.03  0 

0.310 

-0.00807691 

-1.05504028 

0.00205924 

1.05504228 

0.070 

0.310 

“C. 01615686 

-1.04734273 

0.00555283 

1.04735745 

0.110 

0.310 

-0.02116177 

-1.03961727 

0.01037052 

1.03966899 

0.150 

0.310 

-0.02308992 

-1.03182937 

0.01748594 

1.03197751 

0.190 

0.310 

-0.02193956 

-1.02387067 

0.02910578 

1.02428428 

0.230 

0.310 

-0.01770895 

-1.01528662 

0.05152088 

1.01659299 

0.270 

0,310 

-0.01039633 

-1.00262343 

0.11254379 

1.00892013 

0.307 

0.313 

-0.00196843 

-0.90949464 

0.42216929 

1 .00279006 

0.030 

0.270 

-0.00692050 

-1.04718486 

0.00245159 

1.04718773 

0.070 

0.270 

-0.01345874 

-1.03948153 

0 . 00679358 

1.03950372 

0.110 

0.270 

-0.01692233 

-1.03173263 

0.01319399 

1.03181698 

0.150 

0.270 

-0.01730977 

-1.02385556 

0.02364408 

1.02412852 

0.190 

0.270 

-0.01461953 

-1.01549485 

0.0438541 6 

1.01644188 

0.230 

0.270 

-0.00885010 

-1 .00385655 

0.09943140 

1 .00876884 

0.267 

0.273 

-0.00171047 

-0.91949526 

0.39948646 

1.00252729 

0.030 

0.230 

-0.00576523 

-1.03933267 

0.00298813 

1.03933696 

0.07  0 

0.230 

-0.01076326 

-1.03161792 

0.00861383 

1.03165388 

0.110 

0.230 

-0.01268705 

-1.02381455 

0.01779183 

1.02396912 

0.150 

0.230 

-0.01153528 

“1.01566193 

0.03558312 

1.01628506 

0.190 

0.230 

-0.00730669 

-1.00502172 

0.08502685 

1 .00861202 

0.227 

0.233 

-0.00145301 

-0.93096692 

0.37148217 

1.00234644 

0.030 

0.190 

-0.00461092 

-1.03148251 

0.00377684 

1.03148941 

0.070 

0.190 

-0.00807004 

-1.02374182 

0.01158694 

1.02380739 

0.110 

0.190 

-0.00845530 

-1.01577325 

0.02675307 

1.01612548 

0.150 

0.190 

-0.00576564 

-1.00606358 

0.06935426 

1.00845125 

0.187 

0.193 

“0.00119598 

-0.94391469 

0.33667415 

1.00215987 

0.030 

0.150 

-0.00345741 

-1.02363201 

0.00506896 

1.02364456 

0.070 

0.150 

-0.00537869 

-1.01581479 

0.01741257 

1.01596401 

0.110 

0.150 

-0.00422649 

“1.00692324 

0.05244603 

1.00828814 

0.147 

0.153 

-0.00093930 

-0.95816901 

0.29301283 

1.00197024 

0.030 

0.110 

-0.00230454 

-1.01577337 

0.00761441 

1.01580189 

0.070 

0.110 

-0.00268880 

“1.00753898 

0. 03434670 

1.00812424 

0.107 

0.113 

-0.00068291 

-0.97318004 

0.23766570 

1.00178061 

0.030 

0.070 

-0.00115212 

-1.00784771 

0.01511684 

1.00796106 

0.067 

0.C73 

-0.00042672 

-0.98761289 

0. 166771 53 

1.00159470 

0.027 

a. 033 

-0. 00017067 

-0.99858245 

0.07529678 

1.00141724 
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COMPUTER  PROGRAM  C 

DIMENSION  P  (  5 1  *51)  *  T ( 5 1 »  51  )  *TDX(51»51)  *TDY(51*51) 

READ  (5*100)  E*H,N 

NN=N-1 

HH=H*H 

H2  =HH/ 2 • 

H6=H/6. 

DO  5  1=1, N 

5  P( I ,N)=0.0 
Y=-H 

DO  7  J=1 , N 
Y=  Y+H 
X  =  -H 

DO  7  I=1*J 
X=  X+H 
SG=-1 • 0 
WX  =  0  • 

WY  =  0  • 

wxx^o. 

WXY^O. 

DO  6  N!=l ,9 
EN=NI-1 

GG=(FN+0.5)*3.14159265 

GGY*GG*Y 

GGX=GG*X 

SHYssO.  5*(  EXP  ( GGY )  -EXP  ( -GGY )  ) 

CH Y=EXP ( -GGY 1 +SHY 

ANB=0,5*(EXP(GG)+EXP (-GG) )*( l.+2.*FN ) 

ANC=ANB*( l.+2.*FN) 

SG=-1.*$G 
CS=CO$ ( GGX ) 

SN  =  $  IN ( GGX ) 

WX=WX+SG*SHY*CS/ANC 
WY=WY+SG*CHY*SN/ANC 
WXX=:WXX+SG*SHY*SN/ ANB 

6  WXY=WXY+$G*CHY*CS/AN8 
CC=16./3. 14159265/3.14159265 
CB=8*/3. 14159265 
EX=-CC*WX 
EY=2.*X-CC*WY 

EXX=CB*WXX 

EYY=-EXX 

EXY=-CB*WXY 

SI=SQRT(EX*EX+EY*EY) 

T (  I  * J ) = 1 • +  E*S I 

TD  X ( I,J)=E*(EX*EXX+EY*(EXY+2. )  )/SI 

7  TDY( I,J)=E*(EX*EXY+£Y*EYY)/SI 
Y=  1 .0 

J=N  +  1 
10  J  =  J-1 
K  =  d-1 
X=  — H 

DO  20  1=1, K 
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X=  X  +  H 

TA  =  S  IN ( P (  I,J)  )/COS(P( I  * J ) ) 

20  P (  I > J~ 1 )=P(  I ,  J ) +TA* ( R  ( 1  +  1 , J ) ”P ( I *J)  ) ~H*( TDX (  I  ,  J ) ~T A*TDY ( I * J ) ) /T ( I , 
1 J) 

Y=  Y-H 

IF  (J.GT.2)  GO  TO  10 
T  Q  =  HH*H/3  2 • 

Y=  0  • 

DO  45  J  =  2  » NN 
Y=  Y+H 
X  =  0. 

1  =  1 

SX  =  T ( I »  J ) *COS ( P (  I  ,  J)  ) 

5 Y  =  T ( I *  J ) *S I N ( P (  I ,J)  ) 

TQ=TQ+(SX*Y+SY*H/4. ) *H2 
IF  (J.E0.2)  GO  TO  42 
K=  J-l 

DO  40  1=2, K 
X=  X  +  H 

S X  =  T (  I  ,J)*COS(P( I , J)  ) 

S Y  =  T ( I  * J)*SIN( P( I , J)  ) 

40  TQ=TQ+ ( SX*Y+SY*X ) *HH 
42  I  =  J 
X=X  +  H 

S  X  =  T ( I , J)*COS(P( I  *  J )  ) 

S Y=  T ( I  * J)*$IN(P{ I , J)  ) 

SXU=SX*F. /6.+T ( I ,J+1 )*COS<P( I*J+1> )/6. 

45  TQ=TQ+(SXU*( Y+H6)+SY*(X~H6) )*H2 
J=N 
Y=  1  •  0 

YY=l.-H/4. 

X=-H 

DO  50  I  =  1 >  NN 
X=  X+H 
SX=T ( I , J) 

S  Y  =  0  • 

A  A  =H2 

IF  (I.FQ.l)  AA=H2/2. 

50  TQ=TQ+SX*AA*YY 
TQ=TQ+HH*YY/8, 

TQQ=  T  0*8 • 

WRITE  (6*103)  TO© 

60  WRITE  (6*104)  TQ 
80  STOP 

100  FORMAT  (F10.3.FR.8 *  16 ) 

103  FORMAT  <3lHK  DIMENSIONLESS  TORQUE  =*F12.8) 

104  FORMAT  ( 3 l H  ONE-EIGHTH  OF  IT  =,F12.8) 

END 

DATA  CARD  - 

0.100  0.020  51 
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